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ABSTRACT 

We have mapped over 50 massive, dense clumps with four dense gas tracers: HCN 
J = 1 — and 3 — 2; and CS J = 2 — 1 and 7 — 6 transitions. Spectral lines of optically 
thin H^^CN 3-2 and C^^S 5-4 were also obtained towards the map centers. These 
maps usually demonstrate single well-peaked distributions at our resolution, even with 
higher J transitions. The size, virial mass, surface density, and mean volume density 
within a well-defined angular size (FWHM) were calculated from the contour maps for 
each transition. We found that transitions with higher effective density usually trace 
the more compact, inner part of the clumps but have larger linewidths, leading to an 
inverse linewidth-size relation using different tracers. The mean surface densities are 
0.29, 0.33, 0.78, 1.09 g cm-^ within FWHM contours of CS 2-1, HCN 1-0, HCN 3- 
2 and CS 7-6, respectively. We find no correlation of Lm with surface density and 
a possible inverse correlation with mean volume density, contrary to some theoretical 
expectations. Molecular line luminosities L'^^^ were derived for each transition. We 
see no evidence in the data for the relation between L'^^^ and mean density posited 
by modelers. The correlation between L'^^^ and the virial mass is roughly linear for 
each dense gas tracer. No obvious correlation was found between the line luminosity 
ratio and infrared luminosity, bolometric temperature, or the Lj^/Myir ratio. A nearly 
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linear correlation was found between the infrared luminosity and the line luminosity of 
all dense gas tracers for these massive, dense clumps, with a lower cutoff in luminosity 
at Liji = 10^'^ Lq. The L//j-L^^jy^_Q correlation agrees well with the one found in 
galaxies. These correlations indicate a constant star formation rate per unit mass from 
the scale of dense clumps to that of distant galaxies when the mass is measured for dense 
gas. These results support the suggestion that starburst galaxies may be understood as 
having a large fraction of gas in dense clumps. 



1. Introduction 

A full understanding of how stars form, both in the Milky Way and in other galaxies, requires 
an improved picture of the formation of massive stars. Most of the stars in our Galaxy form in 
regions of clustered star formation, and nearly all massive stars appear to form in dense clusters 
(Elmegreen 1985, Carpenter 2000, Lada &: Lada 2003). Massive forming stars provide the only 
obvious signposts of star formation in distant regions of our Galaxy and certainly in other galaxies. 
The star formation prescription most widely used in simulations of galaxy formation and evolution, 
the Kennicutt-Schmidt law (Schmidt 1959, Kennicutt 1998) is based on tracers like Ha (e.g., 
Calzetti 2008), which trace the formation only of relatively massive stars. The signs of formation 
of lower mass stars are mostly invisible at the distances of other galaxies. For all these reasons, a 
better understanding of the formation of massive stars is important. 

However, our knowledge of the massive star formation process is still limited, in both obser- 
vational and theoretical aspects, in contrast to the progress made in the study of low-mass star 
formation. Observational progress has been limited by distance, complexity, and the rapid evolution 
of young massive stars. Most high-mass star-forming regions are very distant from us, normally at 
least an order of magnitude farther than the nearby low-mass star forming regions. When they are 
close enough to study in detail, such as the Orion molecular cloud, we find an extremely complex 
environment marked by high turbulence, high density and opacity (e.g., O'Dell 2005), and multiple 
large-scale motions, including possible explosive events. The fast evolution of the new-born mas- 
sive stars leads to a harsh radiative environment. Together with the jets, winds, and outflows, the 
radiation field rapidly destroys the evidence that we need. 

The confusing observational picture allows quite a variety of theoretical models to exist. Re- 
cently, the models have tended to fall into one of two camps. Bonnell et al. (1997, 1998) and Bonnell 
&; Bate (2005) have begun with a clump containing many Jeans masses, which fragments rapidly 
into small stars that compete to accrete the remaining gas. The other camp, represented by McKee 
& Tan (2002, 2003) and Tan, Krumholz & McKee (2006), argue that such a clump will be com- 
posed of a spectrum of turbulent cores, essentially scaled-up versions of the well-studied cores that 
form low-mass stars. The higher mass stars then form because turbulent pressure produces a much 
higher accretion rate. For an interesting discussion of these different approaches, see Krumholz 
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&; Bonnell (2007). Recent simulations including radiative feedback from forming protostars found 

that the extreme fragmentation found in earlier, isothermal simulations was greatly suppressed, 
allowing more massive stars to form (Krumholz et al. 2007, Urban et al. 2010, Krumholz et al. 
2010). 

In this situation, one positive step to study massive star formation is to make statistical studies 
based on systematic surveys towards a large and well-characterized sample of massive star-forming 
regions. Using suitable tracers, which can reveal the physical and dynamical properties of the 
regions, we can provide the basic constraints that successful theories must satisfy. 

Star formation is not uniformly distributed over molecular clouds, but is restricted to regions 
within molecular clouds where the density exceeds some threshold, roughly n > 10^ or 10^ cm~^ 
(e.g., Evans, 2008). There have been many detailed studies of individual regions of massive star 
formation and surveys of likely regions in tracers of dense gas or dust (e.g., Plume et al. 1992, 
1997; Zinchenko, Mattila, &: Toriseva, 1995; Anglada ct al. 1996; Sridharan et al. 2002; Bcuther et 
al. 2002; Wu &l Evans 2003). These surveys have mostly used signposts of massive star formation, 
such as masers, HII regions, or infrared sources with particular colors. Less biased surveys, using 
infrared dark clouds or blind surveys for millimeter wavelength dust continuum, are beginning to 
become available (e.g.. Price et al. 2001, Williams et al. 2007, SchuUer et al. 2009, Aguirre et al. 
2009, Rosolowsky et al. 2009). These will demand follow-up studies, and the question arises of 
what can be learned from different tracers. 

This work will focus on a large sample of massive star forming regions that are in an early 
phase of star formation, as indicated originally by the association with water masers. We have in 
the past referred to these regions as "cores" , but McKee and co-workers prefer to restrict the use of 
"core" to a smaller region that will form a single or small star group, and they refer to regions that 
will form clusters as "clumps." The expectation is that star-forming clumps will contain multiple 
cores. We adopt the McKee terminology here and refer to these objects as clumps. We emphasize 
however, that their properties are very different from the objects usually going by this name. As 
described by Williams et al. (2000), clumps are traced by extinction or C^^O peaks and have 
modest mean densities of n ~ 10^. Most may not be gravitationally bound. The clumps we are 
studying here are very dense ((n) ~ 10^, e.g. Plume et al. 1997), on average 100 times denser than 
the "cores" in low-mass star forming regions (e.g. Mueller et al. 2002). To keep this distinction 
in mind, we refer to them as "dense clumps." Since most also have high masses, we often refer to 
them as "massive dense clumps." 

The sample has been well studied by multiple CS line surveys (Plume et al. 1992, 1997). A 

subsample of these sources has been mapped in the CS 5-4 transition (Shirley et al. 2003) and 
in 350 lira dust continuum emission (Mueller ct al. 2002). Based on the dust continuum data, a 
ID model was constructed to determine the mass, bolometric luminosity, and density distribution 
of those massive dense clumps; a power-law density distribution (n(r) oc with {p) = 1.8) was 
found to fit the data well (Mueller et al. 2002). The CS 5-4 maps have provided further constraints 
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on the size and virial mass of the dense clumps (Shirley et al. 2003). 

While the maps of CS 5-4 emission and dust continuum emission provided valuable statistical 
data on sizes and masses, systematic studies of other commonly used tracers can provide a valuable 
comparison set for surveys in these tracers, both in distant regions of our Galaxy and in other 
galaxies. In particular, HON has been used much more widely in surveys of starburst galaxies, and 
we need a set of fiducial objects in order to understand those surveys. 

Additional molecular lines that can trace dense clumps, as well as their maps, are needed 
for this purpose. We selected the CS 7-6, CS 2-1, HCN 3-2, and HCN 1-0 transitions for further 
mapping. CS 7-6 and CS 2-1 maps can constrain models of the clumps that have been studied 
with CS 5-4 at different excitation levels. HCN 3-2 is a possible indicator of infall in massive star 
forming regions (Wu & Evans 2003), and it may be able to constrain mass inflow rates needed in 
theories of massive star formation. 

Recent work (Gao & Solomon 2004a, b) shows that the luminosity of HCN 1-0 has a tight and 
linear correlation with infrared luminosity, for both nearby normal galaxies and distant starburst 
galaxies, suggesting that star formation efficiency may be constant in these different systems when 
measured against tracers of the dense gas. We can directly test this L^(^^-L//j correlation in 
Galactic dense clumps, by mapping HCN 1-0 towards massive dense clumps. Our results, based on 
the data in this survey, did indeed suggest that the ratio of massive star formation to the massive 
dense gas is similar in Galactic dense clumps as in starburst galaxies (Wu et al. 2005). That result 
suggests that HCN can be used to connect star formation on scales ranging from individual dense 
clumps in our Galaxy to distant galaxies. Then we may be able to apply our knowledge of Galactic 
massive star formation obtained from CS and HCN studies to star formation in other galaxies. 

For all these purposes, and for building a complete dataset to constrain models of individual 
massive dense clumps, we have carried out mapping surveys using HCN 3-2, 1-0, CS 7-6 and 2-1. 
In this paper we report the basic results of the survey. 

2. Observations 
2.1. The Sample 

The sample in the survey is a subset of the larger sample of dense clumps associated with H2O 
masers that have been surveyed with several CS transitions (Plume et al. 1992, 1997), most of 
which have been mapped with CS 5-4 and dust emission (Shirley et al. 2003, Mueller et al. 2002). 
Table [T] lists the information on sources that have been mapped in this survey. The sources mapped 
in CS 5-4 (Shirley et al. 2003) have virial masses within the nominal core radius (Res) ranging 
from 30 Mq to 2750 M©, with a mean of 920 Mq. The sources in this category have infrared 
luminosities ranging from 10^ - 10'' Lq and most contain compact or ultracompact H H (UCHH) 
regions. These dense clumps are mostly quite massive and have been well studied by several tracers; 
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we will focus on these sources for statistical studies of the properties of massive dense clumps. To 
extend the sample towards lower luminosities for the study of L^^^-L/^j correlation, we selected 
14 IRAS sources from outflow surveys (Zhang et al. 2005; Wu et al. 2004) and a few lower mass 
clumps from other publications. These added sources are listed in part 2 of Tabled! separated by 
a line from the well studied massive clumps in Tabled) These clumps have various luminosities; in 
this paper they are not included in our statistics of the properties of massive dense clumps. 

We mapped massive clumps in the HCN 3-2, HCN 1-0, CS 7-6, and CS 2-1 transitions. Because 
of observing constraints, maps of a particular transition are not available for all sources in our 
sample, but we have attempted to map as many sources as possible in all the tracers. In addition 
to mapping the clumps in the main transitions, we also observed the optically thin isotopes of 
HCN and CS (H^^CN 3-2 and C^^S 5-4) at the center position of most of these massive clumps to 
constrain models and to study infall in massive clumps. Observations were taken between 1996 to 
2007, with information presented in Tabled 



2.2. Observations with the FCRAO 

Maps of HCN 1-0 and CS 2-1 were made with the 14-m telescope of the Five College Radio 
Astronomy Observatory (FCRAO). The 16-element focal plane array (SEQUOIA) was used, with 
typical system temperatures 100-200K. A velocity resolution of 0.1 km s^^ was achieved with the 
25 MHz bandwidth on the dual channel correlator (DCC). We converted the measured to Tr via 
Tr = T^/iriFSSiJc), with rjpss = 0.7. The value of t]c depends on source size; for the typical map 
in this study (~ 10'), r]c = 0.7. The map size was extended until the edge of the HCN J = 1 ^ 
and CS J = 2 ^ 1 emission was reached, typically at the 2a level (mean a ~ 0.3 K km s~^), so we 
could obtain the total line luminosity. 



2.3. Observations with the CSO 

All the other observations were made at the 10.4-m telescope of the Caltech Submillimeter 
Observatorj0 (CSO). The observing date, line frequency, beam size, and main beam efficiency of 
each line are listed in Table 2. The rest frequencies of HCN 3 - 2, H^^qj^ 3 _ 2^ and C^^S 5-4 hues 
have been updated (Ahrens et al. 2002, Gottlieb et al. 2003) after some of our observations. We 
have corrected our data to the new frequencies listed in Table 2. The position switching mode was 
used; reference positions were checked when necessary. This is important for tracers that may show 
absorption features in the spectra (for example, HCN 3-2). We checked all the reference positions 
for HCN 3-2 emission and found them to be free from emission (T^ < 0.5 K). Pointing was checked 



^The CSO is operated by the California Institute of Technology under funding from the National Science Foun- 
dation, contract AST 90-15755. 



- 6 - 



periodically using planets and CO-bright stars. The pointing accuracy was usually better than 6" 
for all the runs. Map sizes of HCN 3-2 and CS 7-6 were extended to at least the point where the 
emission was one third of the peak strength, but usually maps extend to as low as one tenth of 
the peak strength (typically the 3a level, with the mean u ~ 1 K km s"^). The beam sizes and 
efficiencies were fairly constant for all transitions except for the CS 7-6, for which various different 
receiver optics were in place at different times. All data were scaled by the efficiency for the run 
on which they were taken to obtain the values of given in tables and used in the analysis. 

3. Observational Results 

Maps were obtained for 56 clumps with HCN 1-0, 43 clumps with HCN 3-2, 56 clumps with 
CS 2-1, and 51 clumps with CS 7-6. The spectral line information at the central (0,0) position of 
each map is listed in Table [3] to Table [H 

3.1. CS Results 

Table [3] presents the peak temperature of the line (TJj), integrated intensity (jT|jdi^), LSR 
velocity, and FWHM hnewidth of CS 2-1, CS 7-6, and C^^S 5-4 transitions for massive dense 
clumps. Central position spectra of CS 2-1 and CS 7-6 have also been observed by Plume et al. 
(1992, 1997), who obtained CS 2-1 lines with the IRAM 30-m telescope at Pico Veleta, and CS 7-6 
lines from the CSO. Comparing to the results in Plume et al. (1992, 1997), the observed parameters 
of CS 7-6 are consistent for the two CSO observations. The mean ratio of FWHM linewidth in the 
current observations to those from Plume's earlier data is l.libO.4, with median ratio of 1.1. The 
mean ratio for jT'^du is l.libO.7, with a median ratio of 1.0. For the results of CS 2-1 observations, 
the IRAM data give a stronger T|j and larger jTJjdi/; the ratio of jTJjdi^ for the IRAM data to 
current FCRAO data is 3.3ibl.8. This difference comes from the different telescope sizes and beam 
sizes of the IRAM and the FCRAO. But the linewidth from the two observations are close: the 
FWHM ratio of IRAM CS 2-1 data to FCRAO data is 1.1±0.2, and the median is 1.1. 

3.2. HCN Results 

Table m presents the line information for HCN 3-2 and its isotopologue H^^CN 3-2. The mean 
integrated intensity ratio of HCN 3-2 to H^^CN 3-2 is 9.7ib6.8. This ratio is smaller than the 
usual ^^C/^^C isotopic abundance ratio, partly due to the large fraction of self-absorption features 
presented in HCN 3-2 lines. About half of the sources in the massive core sample show double 
peaks. The description of line profiles is in Table [H 

Table [5] lists the observational results for HCN 1-0. HCN 1-0 has three hyperfine components 
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(F=l-0, F=2-l, F=l-1), separated by -7.064 km s"^ (F=l-0 to F=2-l) and 4.842 km s-^(F=l-l to 
F=2-l). We tried to fit the spectra with 3 gaussian components with the same hne width. The peak 
temperatures of every component, the LSR velocity of the main component (F=2-l), the hnewidth, 
and the integrated intensity are presented in Table [5l Massive clumps usually have large linewidths; 
sometimes the line width is wider than the separation of hyperfine components, so that two or three 
of these components will overlap. Sometimes one or two hyperfine components are optically thick 
enough to show absorption. We have noted these features in Table O If the HCN 1-0 line is 
optically thin and the hyperfine levels are populated according to LTE, the line ratios between the 
three hyperfine components are I(F=1-0):I(F=2-1):I(F=1-1)=1:5:3. The observed ratios for most 
of our clumps (Table [5]) do not agree with these optically thin ratios, as expected for lines that are 
almost surely optically thick. This hyperfine anomaly has also been reported in other massive star 
forming regions (e.g., Pirogov et al. 1996, Pirogov 1999). 



3.3. Morphology and Multiplicity 

Contour maps of HCN 1-0, HCN 3-2, CS 2-1, and CS 7-6 are presented in Figures [l]l23 These 
maps show a variety of morphologies, with some multiple peaks. HCN 1-0 and CS 2-1 maps cover 
a much larger area (typically greater than 7' x 7') than the higher J transition maps. A substantial 
fraction of these low-J transition maps show multiple peaks. For HCN 1-0, 16 out of 53 maps (30%) 
show multiple peaks; For CS 2-1, 21 out of 53 maps show multiple peaks (40%). This ratio is lower 
for high-J HCN and CS tracers. Only 8 out of 42 HCN 3-2 maps (19%) and 1 of 51 CS 7-6 (2%) 
present multiple peaks. For CS 5-4 maps (Shirley et al. 2003), this ratio is 16% (9 out of 55). The 
majority of HCN 3-2 and CS 7-6 maps show a single, well-peaked distribution of intensity, similar 
to earlier maps of the CS 5-4 transition or the 350 fim dust continuum emission (Shirley at al. 
2003, Mueller et al. 2002). 

Considering only sources with maps of multiple lines of the same species, there are 37 dense 
clumps mapped in both HCN 1-0 and HCN 3-2; 6 clumps that are single clumps in HCN 1-0 maps 
break into multiple clumps in HCN 3-2. Among the 38 clumps mapped by both CS 2-1 and CS 
5-4, 4 single clumps in CS 2-1 maps show multiple components in CS 5-4 maps. This effect may 
be due to the higher resolution of the HCN 3-2 and CS 5-4 observations. For the 4 clumps that 
have been resolved by CS 5-4 but not CS 2-1 maps, they are either not observed by CS 7-6, or the 
other components are beyond the CS 7-6 mapping area. Consequently, we do not know if they are 
also separate peaks in the CS 7-6 maps. 

The overall impression is that low excitation lines may show secondary peaks that are not 
dense enough to show up in the lines requiring higher excitation, while maps with higher resolution 
may find multiple peaks within the low-resolution map. Because lines of higher excitation in our 
study have better spatial resolution, those multiple peaks may not be seen in the lower excitation 
tracers, just because of instrumental limitations. To check this possibility, we have smoothed the 
HCN 3-2 maps for four clumps with multiple peaks to the resolution of the HCN 1-0 data and 
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compare them to the HCN 1-0 data (Figure The similarity of the smoothed HCN 3-2 maps to 
the HCN 1-0 maps supports the suggestion that HCN 1-0 maps with higher resolution would also 
resolve multiple peaks. 

These few clumps aside, the general appearance of most maps at the current resolution is 
not one of many cores within a clump, but of single well-peaked distributions. Maps with higher 
resolution that still retain sensitivity to extended structure are needed to study substructure in 
these clumps. 



3.4. Clump Sizes 

The size of the clump is characterized by the angular diameter or linear radius after beam 
deconvolution. OtransiUon is the angular diameter of a circle that has the same area as the half peak 
intensity contour: 

n — 9(^^111 — ^ beam \l/2 

"transition — ^\ h / ' 

vr 4 

and Rtransition = GtransitionD /2, where A1/2 IS the area within the contour of half peak intensity, 
(^beam IS the FWHM beam size, and D is the distance to the source. It is important to note that 
these sizes do not represent sharp boundaries. To show this fact, we chose GlO.6-0.4, W3(0H) and 
S231 as examples, since these sources are round in the contour maps, thus better approximated by 
a ID radial profile. In Figure the emission generally declines smoothly to the noise level, when 
plotted against the angular source size, in units of the FWHM angular extent. 

The calculated sizes of clumps in each transition are listed in Table El with the mean and 
median values given at the end of the Table. The statistics on median and mean value include 
only the massive sources, but not those sources below the horizontal line in the tables, which are 
mostly not very massive and were added to study the lower luminosity end of the star formation 
law. The histogram of the deconvolved FWHM sizes (angular size and linear size) of each transition 
are presented in Figures [26] to [271 iii which again we only include massive dense clumps above the 
horizontal line in the tables. Figure [26] shows the histogram of deconvolved angular FWHM size 
of these dense clumps, with a dashed line indicating the beam size in that transition. In HCN 3-2 
maps, 29% of dense clumps have their deconvolved FWHM sizes smaller than the beam size. For 
CS 7-6 maps, this ratio increases to 55%, including one source (G12.21-0.10) not being resolved at 
all. The beam size effect can be nearly ignored for lower-J transitions; almost all the FWHM sizes 
of HCN 1-0 and CS 2-1 are well above their beam sizes. 

As Figure [271 makes particularly clear, the size distribution is skewed to low values before it is 
cut off by our resolution limits. Note that the plots of the lower excitation lines use a larger scale 
because they cover a substantially larger range of sizes. 
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3.5. Luminosities 

For later comparison to work on other galaxies, we compute the line luminosities and the 
infrared luminosities in the same way as is done for extragalactic studies. The line luminosity of 
each core, assuming a gaussian brightness distribution for the source and a gaussian beam, is (we 
use HCN 1-0 for example) 

L'hcni-o = 23.5 X 10-6 X x i^^) x ( ^'^f— ) x J TndvKkm s'K (1) 

Here D is the distance in kpc, and and Obeam are the angular size of the half-peak HCN 1-0 
contour and the beam in arcsecond. This method is analogous to that of Gao & Solomon (2004b), 
but adapted to Galactic clumps. 

The total infrared luminosity (8-1000 fim) was calculated based on the 4 IRAS bands (Sanders 
and Mirabel 1996), as was done for the galaxy sample of Gao & Solomon (2004a): 

Lin = 0.56 xD^x (13.48 x /12 + 5.16 x /as + 2.58 x /go + /loo), (2) 

where fx is the flux in band x from the four IRAS bands in units of Jy, D is in kpc, and L//j is in 
L0. Some authors used the bolometric luminosity (L^o;), which is calculated from the SED of the 
source (e.g. Mueller et al. 2002). Our statistics show that L^o; calculated in Mueller et al. (2002) 
is very close to the L/^j calculated in this paper for the same clumps. The ratio of L/^j to L^oi is 
1.28 ±0.90, and the median is 1.09. The relatively large dispersion for this ratio is strongly affected 
by two sources (S252A and CepA) whose ratio is much larger than others. Excluding these two 
sources, the mean ratio becomes 1.10±0.25, with the median 1.08. The mean and median values 
of Liji are (4.7± 1.2) x 10^ L© and 1.06 x 10^ L©, with values ranging from 708 to 5 x 10^ L©, 
considering only the sources in the original sample. 

The infrared and line luminosities are tabulated in Table [3 with mean and median values for 
the massive, dense clumps at the bottom. The histograms of luminosity are shown in Figure [28l 
While the histograms appear to be well-peaked, we are systematically biased against low luminosity 
sources. Together with the fact that the histograms plot the logarithm of the luminosity, the highest 
luminosity sources are clearly rare. For example, the mean value of luminosity is typically 3 times 
the median. 

4. Discussion 

4.1. Trends in Size with Tracer 

Although all these lines trace "dense" gas, there are differences. Higher J transitions of a 
given species trace denser parts of the core. Crudely speaking, for a given J level, HCN transitions 
have higher critical densities than do the corresponding transitions in CS. The dependence on 
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temperature also differs, with the upper level of the CS 7-6 transition lying 66 K above ground, 
compared to only 26 K above ground for the HCN 3-2 transition. 

In Table [8l we list the properties of the transitions and the statistics of the size ratios and 
luminosity ratios of maps derived from different tracers. Table [9] shows that maps of low-J transi- 
tions have larger deconvolved sizes than those of higher J transitions. For instance, the linear size 
of the HCN 3-2 map is 0.28 of the size of the HCN 1-0 maps, and the ratio is the same for CS 7-6 
versus CS 2-1. 

This result would not be found in a "raisin-pudding" model of clumps: a roughly uniform 
pudding (the clump) with multiple dense cores sprinkled randomly throughout. That model would 
predict similar sizes for maps of transitions sensitive to different density regimes. Such a pattern 
was seen in early studies of multiple transitions in a small number of sources (Snell et al. 1984, 
Mundy et al. 1986, 1987). The current results have higher spatial resolution and more complete 
maps in the high-excitation lines than the early results and do show that higher excitation lines 
are more centrally condensed. If the clumps are centrally condensed, with either a smooth increase 
in density with decreasing radius, or an increased fraction of the volume filled by dense "raisins", 
the differences in size would be predicted as the lines from higher-J transitions would be strongly 
enhanced in the denser regions. Detailed models of individual dense clumps do support density 
gradients with power laws providing good fits (n(r) oc with {p) = 1.8) (Mueller et al. 2002). 

4.2. An Inverse Size-Linewidth Relation 

The linewidth at the peak of the dense clumps for different tracers shows an interesting trend. 
As presented in Table [TOl the mean and median FWHM linewidth of the dense clumps increase for 
higher J transitions, from about 4.8 km s~^ for CS 2-1, to about 5.9 km s~^ for CS 7-6. If we look 
at a optically thin tracer, like C^^S 5-4, the mean FWHM is 4.7 km s~^, which is closer to that of 
the CS 2-1 than higher J CS transitions. 

Studies of molecular clouds have established a positive correlation of cloud size, as mapped with 
a single tracer, and linewidth of that tracer (Larson 1981, Myers 1985), which can be interpreted 
in terms of the virial theorem. The linewidth-size relation within individual cores or clumps is 
much less well-established (e.g., Evans 1999). Here we refer to the useful distinction introduced 
by Goodman et al. (1998) between 4 different kinds of relationships. A Goodman Type-4 relation 
plots linewidths of a single tracer measured over regions of increasing radius; such plots show little 
correlation (Goodman et al. 1998, Barranco &: Goodman 1998). Goodman Type-1 relations plot 
sizes measured in multiple tracers for multiple cores or clumps. Such plots generally show positive 
correlations when using tracers of modest density applied to low-mass and intermediate-mass cores 
or clumps; for example Caselli &: Myers (1995) found /S.v oc r"'^^^"'^ with a correlation coefficient 
of 0.81 using NH3, CS 2-1, and C^^O 1-0 maps of 8 starless cores. 

In contrast, studies of the cores in this sample with Goodman Type-2 relations (single tracer. 
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multicore) have shown no correlation between Hnewidth and size (Plume et al. 1997 for CS 7-6) 
or a weak (r = 0.43) correlation (Shirley et al. 2003 for C^^S 5-4). In both cases, the linewidths 
were 4 to 5 times larger than would be predicted from the usual linewidth-size relations developed 
from studies of clouds and low mass cores. With maps of multiple tracers toward many cores, we 
could construct Goodman Type 3 relations (multitracer, single core) for each core. Instead, we try 
a modified version (Type 5?) in which we use the mean value over many cores for each transition. 

Figure [29] shows the mean linewidth versus the mean size of the dense clumps for CS 2-1, CS 
5-4, CS 7-6 and HCN 1-0 transitions. HCN 3-2 lines have serious self-absorption, therefore they 
are excluded from this figure. HCN 1-0 has 3 hyperfine components; we took the mean linewidth 
of the strongest component (F=2-l) for the plot. The mean linewidth of the tracer has an inverse 
correlation with the mean size of the dense clumps. For example, while the CS 2-1 sizes are on 
average 4 times larger than those of CS 7-6, the mean of the linewidth ratio of CS 7-6 to CS 2-1 
is 1.18±0.34 for the massive dense clumps; the median ratio is 1.16. The CS 2-1 emission traces 
a larger, quieter area than higher J transitions, while CS 5-4 and CS 7-6 trace more compact and 
apparently more turbulent regions. Alternatively, the larger linewidths toward the center could in 
part reflect higher inflow motions. A similar trend was seen in observations and models of G10.6 
(Keto 1990). It is also possible that the turbulence (here refers to all unresolved gas motions, 
including inflow, outflow, complex fluid motions...) has radial structure, which is higher close to 
the forming protostars and therefore is being picked up better by the higher- J, dense gas tracers. 
The usual size-linewidth relations certainly do not apply to these regions where massive stars are 
forming. 



4.3. What Physical Conditions Do the Lines "Trace"? 

While we speak of a transition as "tracing" gas of a certain density, this formulation is too 
simplistic; the emission from a transition increases smoothly and strongly over a wide range of 
densities and then levels out just above the critical density. [This description applies only to lines 
with frequencies high enough that radiative transitions stimulated by cosmic background radiation 
or continuum emission from dust are unimportant (Evans 1989).] Nonetheless, it is common to 
speak of a transition as tracing its critical density. If anything, the opposite is true, as transitions 
fail to provide useful probes of density above their critical density. A more sensible fiducial density 
is the effective density (Evans 1999), which is the density needed to produce a line with a radiation 
temperature of 1 K for typical conditions. For these comparisons, we use the critical and effective 
densities calculated for a kinetic temperature of 100 K (Evans 1999). Both these measures are 
given in Table [H 

The critical densities (ucrit) of the HCN 1-0 and CS 2-1 lines differ by an order of magnitude, 
but their effective densities (ug//), are very similar. The similarity of their mean sizes indicates that 
n^ff better reflects the excitation requirements. For HCN 3-2, CS 5-4, and CS 7-6, their sizes are 
also similar while the critical and effective densities make opposite predictions; CS 7-6 has a much 
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higher T^eff than the other two, but HCN 3-2 has the highest Ucrit- The excitation energy of the 
upper level (Eup) is considerably higher for the CS 7-6 transition. If we consider only the sources 
that have been mapped by all three high-J transitions CS 7-6, CS 5-4, and HCN 3-2 (column 5 
and 7 in Table [9]) , statistics show that the mean size and luminosity of the dense clumps increase 
with decreasing E^p of the transitions, as well as decreasing ngjj, while the correlation with n^rit is 
weaker. The tracers that need higher excitation energy or density are confined to smaller regions in 
the dense clumps, consistent with centrally condensed structures. In addition, Ug// better reflects 
the excitation requirement than Ucrit- 

Empirically, we find that the HCN 1-0 and CS 2-1 lines form a matched set, while HCN 3-2 
and CS 5-4 lines make another pair; each pair has very similar sizes and fairly similar luminosities 
in the mean. The CS 7-6 lines are close to the properties of HCN 3-2 and CS 5-4 pair, but still 
more compact in size. 



4.4. Virial Mass, Surface Density, Volume Density, and Relations with Ljr 

The virial mass of the dense clumps with power-law density distributions (n(r) oc r"^) is 
calculated by 

, , bRAv'^ (i?/lpc)(At;/lkm s-i)2 

Mv^r{R) = ^ ^ « 209 ^ ^ ^ ' '-Mq (3) 

Qaia2LrLn2 aia2 

where ai is the correction for a power-law density distribution and 02 is the correction for a 
nonspherical shape (Bertoldi & McKee 1992). For aspect ratios less than 2, 02 ~ 1 and can be 
ignored for our sample. To calculate ai we take the average p value (1-77) from Mueller et al. 
(2002) for all massive dense clumps. A detailed discussion of calculating the virial mass of dense 
clumps in this sample can be found in Shirley et al. (2003). 

The virial masses calculated from each transition {Mvir{Rtransition))-, for HCN 1-0, 3-2, CS 
7-6, and 2-1, are tabulated in Table [TTl Rather than using the linewidth of the tracer, we used the 
linewidth of C^^S 5-4 for all these tracers because it is more reliably optically thin than is the 3-2 
line of H^^CN. For lower mass sources (those below the horizontal line in Table [TT]) . their C^^S 5-4 
lines are too weak to detect or to give a Gaussian fit except for IRAS20188+3928. The FWHM 
linewidth of CS 2-1 is closest to that of C^^S 5-4 in our sample, with a ratio of linewidth for CS 
2-1 to C'^^S 5-4 of 1.2ib0.1 and a median of 1.1. So, for these lower mass clumps, we used the CS 
2-1 linewidth to calculate their virial mass, and assigned an additional 20% uncertainty to their 
linewidth uncertainties. While the linewidths vary among transitions by about 25%, some of this 
variation must be caused by optical depth differences, and the differences in any case, are small. 
The sense of these calculations is then an estimate of the mass within the FWHM contour of each 
tracer. The mean ratio of M^r for CS 2-1 to CS7-6 is 4.0, with a median of 3.7; the mean ratio 
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for HCN 1-0 and 3-2 is 2.9, with a median of 2.8. One would expect these ratios to reflect the size 
ratios, which are 4.1 and 3.0 (the mean ratio), respectively. 

The virial mass enclosed by the FWHM contour of dense gas tracers represents the mass, most 
of which is dense gas, that is available to clustered star formation. Therefore, one should expect a 
positive correlation between the star formation rate, as indicated by the infrared luminosity, and 
the virial mass. One caveat is that the star formation rate can be underestimated by L//j until 
the initial mass function is well sampled by the star formation process (Krumholz & Thompson 
2007) . We plot Mvir versus L//j for four transitions in figure [30l It is obvious that the infrared 
luminosity is positively correlated with the virial mass for these clumps. We used linear fits based 
on least squares minimization and alternatively on robust estimation, in which the absolute value 
of the deviation is minimized. Comparing these two fits usually provides a more realistic appraisal 
of the true uncertainties. 

As noted by Wu et al. (2005), there is a cut-off around L'^ IO^'^Lq in the \jjji-L'jj^qj^^_q 
correlation. If we fit only the data for massive dense clumps with L/r > 10^'^ L©, the fitting 
results are as follows (r is the correlation coefficient): 

HCN 1-0: 

Least Squares : log{Ljji) = 1.29(±0.08) x log{Mvir{RHCNi-o)) + 0.98{±0.28);r = 0.69 
Robust Fit : log{LiR) = 0.70 x \og{Mvir{RHCNi-o)) + 3.03 

HCN 3-2: 

Least Squares : log{LjR) = 1.23(±0.09) x log{Mvir{RHCN3-2)) + 1.98(±0.24); r = 0.61 
Robust Fit : log{LiR) = 0.79 x \og{Mvir{RHcm~2)) + 3.16 

CS 2-1: 

Least Squares : log{LiR) = 1.13(±0.06) x \og{Mvir{RcS2-i)) + 1.61(±0.19); r = 0.67 
Robust Fit : log{LjR) = 0.98 x log{Mvir{RcS2-^i)) + 2.24 

CS 7-6: 

Least Squares : log{LjR) = 1.42(±0.07) x log{Mvir{RcS7~e)) + 1.28(±0.23); r = 0.49 
Robust Fit : log{LjR) = 0.62 x log{Mvir{RcS2-i)) + 3.76 

There are substantial variations between fitting method employed and among the line used. If 
we average the two fitting methods for each line and then average over all lines, the resulting slope 
is 1.02 zb 0.03. Thus, Ljr is roughly linearly correlated with the virial mass for dense clumps with 
LjR > 10^-5 L0. 
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Below this threshold, the slope is steeper. We also provide linear fits and robust fits to the 
sources with L/^j < 10^'^ Lq (the CS 7-6 maps have too few sources below this luminosity to make 
a reasonable fitting): 

HCN 1-0: 

Least Squares : log{LjR) = 2.02(±0.13) x log{Mvir{RHCNi-o)) - 1.43(±0.30); r = 0.87 
Robust Fit : log{LiR) = 1.71 x \og{Mvir{RHCNi~o)) - 0.81 

HCN 3-2: 

Least Squares : log{LjR) = 2.30(±0.25) x log{Mvir{RHCN3-2)) - 0.98(±0.53); r = 0.83 
Robust Fit : log{LiR) = 0.97 x \og{Mvir{RHcm-^2)) + 1-69 

CS 2-1: 

Least Squares : log{LiR) = 2.46(±0.23) x \og{Mvir{RcS2-i)) - 2.44(±0.56); r = 0.87 
Robust Fit : log{LjR) = 1.59 x log{Mvir{RcS2-i)) - 0.68 

When only clumps with Ljr < 10^'^ are included in the fit, the Ljr increases steeply with the 
mass of dense gas with a power law ranging from 1.5 to 2.5. Using the same averaging method 
yields a mean slope of 1.84 it 0.20. 

The surface density of a molecular clump should have an important impact on the time scale 
and mass accretion rate of star formation in the dense clump (McKee & Tan 2002, 2003). Recent 
work by Krumholz & McKee (2008) proposed that only clouds with column densities of at least 1 
g cm~^ can form massive stars. We can test this idea on these dense clumps with different dense 
gas tracers. We calculate the surface density within the FWHM size of each transition, following 
the method used by Shirley et al. (2003) for CS 5-4: 

^ _ MyiriRtransition) _ ^ ^^^ {M^jr/l-O X IO^Mq) 

^transition — ~ U.DDO jj- j- ^2 (."J 

-^transition [^transition/ l^PC) 

The calculated surface densities are listed in Table [TTJ The surface densities are lower than 1 g 
cm~^ within the FWHM size of CS 2-1 and HCN 1-0 maps. For massive dense clumps, the mean is 
0.29 g cm~^ for HCN 1-0, with a median 0.28 g cm~^, and 0.33 g cm^^ for CS 2-1, with a median 
0.31 g cm~^. Almost no source has surface density greater than 1 g cm~^, when measured by HCN 
1-0 or CS 2-1. Higher J transitions trace more compact, denser regions, resulting a larger surface 
densities. The mean surface density within the FWHM contour of HCN 3-2 is 0.78ib0.55 g cm~^, 
with a median of 0.68 g cm^^. This result is similar to the result from CS 5-4 maps (the mean is 
0.82 g cm~^, with a median of 0.60 g cm~^ (Shirley et al. 2003). CS 7-6 traces the most compact 
part of the clumps in this survey, and its surface density is closest to the theoretical threshold of 1 
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g cm~^. The mean surface density within CS 7-6 clumps is 1.09 g cm~^, with a median of 0.89 g 
cm^^. The percentage of sources that have surface density larger than 1 g cm~^ is 26%, 31%, and 
38% for clumps traced by CS 5-4, HCN 3-2, and CS 7-6, respectively. 

As one might have expected, tracers with higher n-eff trace higher surface densities. Do the 
surface densities correlate with other properties? One might expect that the higher the surface 
density, the larger the virial mass of the clumps within the same FWHM area. According to 
Krumholz & McKee's (2008) model, a surface density of 1 g cm~^ allows the gravity of the core 
to overcome the radiation pressure of the new born stars, in order to keep the accretion going and 
form massive stars inside. If the surface density of a dense clump is high, one may also expect that 
the star formation rate indicator L//j and the indicator of star formation efficiency (L/^/Myj^) are 
also high. 

We check for any correlations between the surface density and Myir, Ljr, Li^/Myir, and line 
luminosity L' for the four dense gas tracers (HCNl-0, HCN 3-2, CS 2-1, CS 7-6) in Figure ED to [Ml 
There are no obvious correlations found between the surface density and these parameters for dense 
clumps, except for the very lower end of the luminosity for a few low-mass cores, whose surface 
density is lower than massive clumps. In particular, L//j ranges over four orders of magnitude 
within a total variation of Tihcn2,-2 of about a factor of 10 for dense clumps. We see no evidence 
for a sharp threshold in surface density for massive dense clumps. 

Of course, the threshold is probably not exact, and our values have substantial uncertainty, so 
generally speaking, these massive clumps are in the right range to form massive stars. With density 
gradients, most of the dense clumps in our sample will reach the critical value of 1 g cm~^ at some 
radius. The amount of mass above the threshold might help to determine the mass of the stars 
that can form. The CS 7-6 line comes closest to being a tracer of surface density of 1 g cm~-^. The 
mass inside the CS 7-6 contour is on average about 0.25 of the mass inside the CS 2-1 contour and 
has a mean value of 1300 Mq. Thus, these clumps on average have plenty of mass at high enough 
surface densities to form clusters of massive stars. 

The volume density of a cloud is often taken to control the star formation rate of the cloud. We 
calculated the mean volume density (n = n(H2) + n{He) + ... = p/fimn, with fj, = 2.37 (Kauffmann 
et al. 2008) within the FWHM size of the clumps for each transition, according the the following 
equation: 

^viriRtransitio-n) ^ , ^.q i,-^vir / ■^q) /p\ 

^transition — 4 ^ 77? 7l Ts ^ ' 

'S'^f^'I^H ^transition {transition/ i-PC) 

The derived volume densities for each transition are tabulated in Table [TTJ 

The logarithmic averages over all clumps of the mean volume density (logn) are 4.1ib0.1, 
4.3±0.1, 5.1±0.1, and 5.4±0.1 for HCN 1-0, CS 2-1, HCN 3-2, and CS 7-6 respectively. The 
uncertainties are the standard deviation of the mean. In Plume et al. (1997), the volume density 
of a similar sample of clumps were derived by fitting multiple CS and C'^^S spectral lines towards 
the center position with an LVG model for the excitation. The logarithmic average over the sample 
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was (lognLvc) = 5.9 ± 0.2. 

The excitation densities {tilvg) derived from the LVG model (Plume et al. 1997) are larger 
than mean densities (n) derived from the mass and size for two reasons. The value of n is an 
average over all the gas within the FWHM contour traced by a certain tracer. Even the CS 7-6 
maps cover more than one to a few beam sizes. The density derived from single point spectrum 
with an LVG model was calculated from one beam measurement at the peak. Since the clumps 
are centrally condensed, n will be less than n from the LVG model. In addition, the excitation 
requirements of the hardest-to-excite transition in the LVG fit push the fit to the higher densities 
in the region. The fact that the values of n from the highest excitation lines are closest to the 
values from the LVG models is a reflection of this point. These points were discussed in detail by 
Plume et al. (1997). 

If the star formation rate is proportional to the amount of gas and inversely proportional to 
the free-fall time, a simple prediction is that oc n/tfj oc n^'^. This idea has been used [e.g., 
Krumholz et al. (2007), Narayanan et al. (2008)] to explain the basis of the observed Kennicutt- 
Schmidt relations for star formation in galaxies. Is the star formation rate proportional to the 
mean density to the power of 1.5, as theorists have suggested? We plot the mean volume density 
of each tracer vs. infrared luminosity in Figure [35l Linear least square fits with uncertainty give 
the following results: 

Least Squares : log{LiR) = -3.42(±0.28) x \og{n{RHCNi-o)) + 20.11(±1.26), r = -0.55 

Robust Fit : log{LiR) = -1.02 x \og{n{RHcm-o)) +8.78 
Least Squares : log{LiR) = -5.28(±0.73) x log(n(i2//c'7V3_2)) + 33.35(±3.99), r 

Robust Fit : log{LiR) = -1.2 x \og{n{RHcm-2)) + 10.96 
Least Squares : log{LiR) = -3.46(±0.25) x log{n{Rcs2-i)) + 20.46(±1.14), r -- 

Robust Fit : log{LjR) = -0.90 x log{n{Rcs2-i)) + 8.35 
Least Squares : log{LiR) = -1.18(±0.51) x log(n(i?c57-6)) + 11.78(±0.30), r -- 
Robust Fit : log{LjR) = -0.70 x log{n{Rcs7-6)) + 8.90 

Far from the prediction that n and hjR are positively correlated, in general, hjR decreases with 
increasing n. (There may be a hint of positive correlation for the few low-luminosity cores with 
adequate data.) Interestingly, all the fits from robust estimation have slopes near —1 with mean 
value over all tracers of —0.96 it 0.21. Since most massive clumps are large, they have lower mean 
densities; this fact probably explains the negative correlation if luminosity is just proportional to 
mass above a threshold density. A similar effect is seen in the plot of L^^^ versus n (Figure [36]) . 
This result is also plausible because L^^^ is almost linearly correlated with Ljr (see §5.2). 



= -0.57 
= -0.49 
= -0.37 
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We can also plot L/^ versus ulvg^ which traces the denser part of a clump with a density 
gradient. The plot of Lj^ vs. nivG (Figure [37|) . reveals a very weak positive correlation with L//j. 
The linear least square fit with uncertainty gives 

log{LiR) = 15.7(±10.1) X \og{n{LVG)) - 87.9(±119.3), r = 0.18 

We see no sign of Lir oc v}'^ for these massive dense clumps, regardless of how density is 
defined. 

We also ploted the correlation between the star formation efficiency (indicated by the infrared 
luminosity devided by rivial mass) and the virial mass for different tracers (Figure [38|) . It shows 
no trend between them 

Comparing the correlations between L//j and virial mass, surface density, and volume density 
for dense clumps, we found that only the virial mass has a strong correlation with Lm] the correla- 
tions between Lm and surface density or mean volume density are weak. It suggests that the most 
relevant parameter for the star formation rate in a dense clump is how much gas above a threshold 
density it contains. The surface density or mean density within the clump are not the controlling 
factors. 



4.5. Do Line Luminosity Ratios Correlate vi^ith Other Properties? 

Table [8] gives the mean, median, and standard deviation of the ratios between the various 
lines. The luminosity ratio of high J to low J transitions could reflect the ratio of very dense gas 
to less dense gas in a clump. They could also depend on the infrared luminosity because a higher 
luminosity will result in a greater fraction of the gas being warm enough to emit in the higher 
J lines. We have tested this latter possibility by plotting line luminosity ratios versus infrared 
luminosity L//j. Figure [39] indicates that only the ratio ^ F^^"*^ shows a weak correlation with L//j. 

^CS5-4 

It can be fltted by a linear correlation: 

log(-f^) = 0.18(±0.05) X log{LiR) - 1.11(±0.26). 

The correlation coefficient r is 0.49. Furthermore, the spread in values of the ratio is modest for 
most ratios (Table [8]). The linear fit for other ratios versus Lir show weak or no correlations: 

^' 

log(-f^^) = 0.09(±0.10) X log{LiR) - 1.28(±0.53), r = 0.17 

^CS2-l 

log( f^^"^ ) = -0.01(±0.06) X log{LiR) - 0.54(±0.30), r = -0.03 

^CS2-l 

log{^l£M^) = 0.084(±0.03) X log{LjR) - 0.10(±0.14), r = 0.46 

^HCNl-O 
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Thus, with the possible exception of the CS 7-6 hne, we find httle dependence of the fine ratios on 
luminosity, suggesting that temperature differences are not important for most ratios. In particular, 
the ratio of ig nearly constant for these clumps, with much less scatter than other ratios. 

^HCN3-2 

These ratios are primarily sensitive to the density structure for these massive clumps, and the uni- 
formity of the values reflects the fact that most dense clumps have power law density distributions 
(Mueller et al. 2002, Shirley et al. 2003). 

We can also correlate the line ratios with the bolometric temperature (Tfeo/), which is the 
temperature of a blackbody with the same mean frequency as the actual SED. In low mass cores, 
Tboi is related to the evolutionary status of a core (Chen et al. 1995), with lower T^oi corresponding 
to earlier stages. In these sources, it might also reflect the relative amounts of gas at higher and 
lower temperatures. The plots in Figure HOl show little correlation between line ratios and Tf,o;. The 
linear least square flts to these correlations give: 

log(-f^^) = 0.005(±0.003) X Tboi - 0.67(±0.22), r = 0.32 

u 

log( J^^"^~^ ) = -0.005(±0.004) X Ttoi - 0.50(±0.32).r = -0.22 
ll 

log(-f^^) = -0.005(±0.003) X T^oi - 0.25(±0.24).r = -0.27 

u 

log( = 0.0006(±0.002) X T^oi - 0.65(±0.20).r = 0.05 

This result suggests that the fraction of gas in various density regimes does not depend strongly 
on evolutionary stage, at least for those stages probed by our sample. 

Do the line ratios correlate with L//j/M„jr, which is a star formation efficiency indicator, if 
Li^ is tracing star formation rate? The plots in Figure [4T] show no clear trends. These results 
suggest that all of these tracers are probing gas that is dense enough to be involved in massive star 
formation. The lower J lines simply trace a larger region of somewhat lower density and thus have 
larger map sizes and luminosities. 

Another approach to an evolutionary sequence for massive star-forming regions is based on the 
evolution of HII regions. In this picture, a clump would evolve from having no radio continuum 
(NoRC), to an Ultra-compact HII region stage (UCHII), to a Compact HII region stage, and finally 
to a complete HII region stage. In the paper of Shirley et al. (2003), an extensive literature search 
using the SIMBAD database was performed to find HII regions associated with these dense clumps. 
Using the taxonomy of Kurtz (2002), HII regions are classified as UCHII if the diameter of 2 cm 
radio continuum emission is <0.1 pc, CHII is the diameter is <0.5 pc, and an extended HII region if 
the diameter is greater than 0.5 pc or clearly associated with a classical HII region. The classification 
of the evolutionary states of these dense clumps refers to Table 1 of Shirley et al. (2003). If we 
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classify the massive clumps in the sample according to this scheme, do we see any trends? We plot 
the line ratio of (the ratio of two tracers that form a well matched sample in our survey, 

so we have enough points to make the statistics) versus star formation efficiency in Figure [l2] with 
different symbols for each HII region status. No obvious differences emerge in the plot. The mean, 
median, and standard deviation of this ratio for each evolutionary stage are tabulated in Table [I2j 
We see no obvious trend for the ratio of along the evolutionary sequence, except that 

the mean and median ratio is lower for sources with no radio continuum emission than other three 
categories, indicating either that the dense gas ratio is lower or the tempertures are lower in this 
stage. 



4.6. Dynamical Properties 

The dynamical state of molecular clumps is an important input into star formation theories. 
As has been discussed many times, the linewidths of massive clumps are much wider than expected 
for thermal broadening, and some form of turbulence is clearly important. For example, the char- 
acteristic temperature in these dense clumps is less than about lOOK. At this temperature, the 
sound speed is 0.7 km s"^. The mean FWHM linewidth from an optically thin tracer C^^S 5-4 for 
these massive sources is 4.7 km s~^. This gives a Mach number of 6.7. Thermal broadening can 
only account for a small fraction of this broad linewidth. The turbulence is clearly supersonic and 
seems to increase on smaller scales, as discussed above. 

Evidence for infalling gas can be seen in the form of line profiles skewed to the blue in optically 
thick lines compared to optically thin lines. Evidence for a significant excess of blue skewed profiles 
over red-skewed profiles was seen in low-mass Class (Gregersen et al. 1997) and Class I sources 
(Gregersen et al. 2000). An excess was also found in a sample of low-mass starless cores (Gregersen 
& Evans 2000). A recent study of somewhat more massive starless cores in Orion found no excess 
of blue-skewed profiles (Velusamy et al. 2008). While Velusamy et al. do not give masses of their 
cores, they are drawn from the sample of Li et al. (2007), which has a maximum mass of 46 Mq. 
These studies all refer to what are now called cores, plausible precursors of single stars or small 
groups. The massive clumps in this study are much more massive and are likely to be the formation 
sites of clusters. In such sources, a blue-skewed profile would not reflect infall onto a forming star, 
but inflow toward the central parts of the clump. 

The turbulence in these massive clumps makes detection of inflow motions difficult, but some 
sources have shown line profiles similar to those that are interpreted as evidence of inflow in lower 
mass cores. In particular, we found an excess of "blue" profiles over "red" profiles in initial ob- 
servations of the HCN 3-2 line for this survey (Wu & Evans 2003), and a similar result was found 
using HCO"'" by Fuller, Williams, &: Sridharan (2005). 

Since we have observed more sources with the HCN 3-2 line and an optically thin line (48) 
since the paper (Wu &: Evans 2003) on the statistics of blue profiles of 28 massive clumps, we can 
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reassess this situation. The definitions of "blue" and "red" profiles have been discussed by Wu & 
Evans (2003) and references therein. Here we note only that inflow may produce an asymmetric 
line profile in an optically thick line, skewing the peak toward the lower velocities relative to the 
peak of an optically thin line. This offset, measured in units of the width of an optically thin line, 
is defined by 



Vthick '^thin ^j-^ 

An absolute value 6v > 0.25 is considered significant. In addition, the line profile may be 
double peaked with the blue peak stronger than the red peak. We refer to this signature as the 
blue peak. Since some profiles show the opposite asymmetries, one can characterize a sample by 
the excess, defined as 



^total 

We find that B^^CN 3-2 is aheady optically thick in some sources that show absorption features, 
while C^'^S 5-4 may be a better optically thin tracer for calculating 5v. The mean ratio of the FWHM 
of H^^cN 3 _ 2 to C^^S 5-4 for our sample is 1.27±0.32. So we use C^^S 5-4 as the optically thin 
line to calculate 5v in this paper, and we use H^^CN 3-2 only when the C^'^S 5-4 line is not 
available. The 6v values are listed in Table 4. Using this measure, we identified 21 blue and 15 
red profiles out of 48 sources; this leads to £^ = 0.13 compared to E = 0.29 in Wu & Evans (2003). 
Using the blue peak method, we identified 16 blue and 9 red profiles and E = 0.15, compared to 
0.29 in Wu & Evans (2003). Thus, the inflow signature is not as common in the full sample as it 
was in the initial results. However, there are still more blue asymmetries than red, suggesting that 
mass inflow is occurring in at least some of the sample. 

For comparison. Fuller et al. (2005) found values of E of 0.15 for the HCO^ 1-0 line and 0.19 
for the H2CO JK^lK-^.l = 2i2 lii line in their study of 112 objects. They found no significant 
excess in the 3-2 or 4-3 lines of HCO"*". There are only 3 sources in both the sample of Fuller et 
al. (2005) and our sample: S231, G59.78+0.06 and IRAS20106+3545 (called WFS2, 79 and 93 in 
Fuller's sample). In both samples, S231 was classified as a red profile and G59.78+0.06 as a blue 
profile (infall candidate). We don't have optically thin line data for IRAS20106, but its HCN 3-2 
line show no asymmetry, which agrees with the observations by Fuller et al. (2005). 

We also tabulate in Table [13] evidence for outflows based on the outflow catalog of Wu et 
al. (2005). Noted that this outflow catalog doesn't cover all the sources in our sample. It is a 
collective catalog that provides all the reported outflows by 2004. Therefore the following statistics 
of outflows are based on the current available outflow reports, but not complete surveys. Usually 
outflow and infall signatures are not present in the same source, but there are obvious exceptions. 
Of the 46 sources whose dv are available, outflow has been reported in 16. Of these, 8 have red 
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asymmetry with Sv > 0.20, 6 have blue asymmetry with Sv < —0.20, in which outflows and infall 

may occur at the same time. If we look the rate of the appearance of the outflow or infall signature 
as the signature of "dynamical activity", the rate of the sources that have infall or outflow to the 
total for each evolutionary stages are 55% for regions without HII regions, 66% for UCHII regions, 
64% for CHII regions, and 33% for fully formed HII regions. These numbers suggest that the 
UCHII/CHII stages are more dynamically active than later stages. 

5. Star Formation Laws: Connecting to Extragalactic Star Formation 

A linear L7R-L^fi^^_Q correlation has been reported (Gao &: Solomon 2004a) for nearby 

and distant starburst galaxies. A similar correlation has been found for Galactic clumps when 
studying CS 5-4 and dust emission maps (Shirley 2003, Mueller et al. 2002). However, the study 
of the Kennicutt-Sclimidt law with CO and HI shows a super-linear correlation (Kennicutt 1998, 
Kennicutt et al. 2007). Why are these two correlations different? 

Two distinct explanations have been offered. Gao & Solomon (2004a) and Wu et al. (2005) 
suggested that the line luminosity traces the mass of dense gas; the linear relation between Lir 
and line luminosity then reflects the fact that the star formation rate is linearly proportional to 
the amount of gas above the threshold for star formation. The fact that massive dense clumps in 
our Galaxy lie on the same relationship seen for starburst galaxies indicates that starburst galaxies 
can be understood as large collections of dense clumps, such as we observe in our Galaxy. For a 
worked example of this idea, see Wu et al. (2009) . 

Krumholz & Thompson (2007) and Narayanan et al. (2008) argue that the underlying star 
formation law is superlinear and the linear relations seen for dense gas tracers result from cancella- 
tion from a superlinear dependence of star formation rate on density. As explained by Narayanan 
et al. (2008), the following two relations. 



SFR (X p 




can be combined to give 



SFR oc L' 



I a 
'mol ' 



with a = N/13. In this picture, the linearity for dense gas tracers (a ~ 1) results from the 
coincidence that and /3 have similar values. They predict that values of a will decline as higher 
excitation lines are used because (5 is larger for subthermal lines. While Narayanan et al. (2008) 
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are somewhat vague about the definition of p, it is generally discussed in terms of mean density 
along the line of sight, hence proportional to surface density. 

The current surveys allow us to examine the correlations for dense clumps in our Galaxy for 
five different tracers of dense gas. For tracers with sufficient data on other galaxies, we can compare 
the Galactic data on individual clumps to that for other galaxies. 



5.1. What do L'^^^ Measurements Trace? 

A first question is whether L'^^j^^^q traces the mass of dense gas, as suggested by Gao &: 
Solomon (2004a). We plot the L'jj(jj^^_q-Mvir{Rhcni-o) correlation in figure 1331 For comparison, 
we also plot the Lmoiecuie vs. My I ii{Rmoiecuie) for the other four dense gas tracers in figure HS) For 
all the panels, we show fits based on least squares minimization (dashed line) and robust estimation 
(solid line). Here are the results from the fits: 

HCN 1-0: 

Robust Fit : logiL'^cm^o) = 1-04 x log{Mvir{RHCNi-o)) - 1-35 

Least Squares : logiL'^cm-o) = 1.17(±0.05) x log{Mvir{RHcm-o)) - 1.78(±0.16); r = 0.91 
HCN 3-2: 

Robust Fit : log{L'HCN3-2) = 1-17 X log{Mvir{RHCN3-2)) - 1.89 

Least Squares : log{L'fjcN3-2) = 1.33(±0.06) x log{Mvir{RHCN3-2)) - 2.36(±0.18); r = 0.93 
CS 2-1: 

Robust Fit : log{L'cs2^i) = 1-06 x \og{Mvir{RcS2~i)) - 1-50 
Least Squares : log{L'cs2-i) = 1.26(±0.05) x \og{Mvir{RcS2-i)) - 2.17(±0.16); r = 0.91 
CS 5-4: 

Robust Fit : log{L'cs->-i) = 0.95 x \og{Mvir{Rcsb~A)) - 1-26 
Least Squares : log{L'c s-,-a) = 1.26(±0.06) x \og{Mvir{RcS5~i)) - 2.23(±0.17); r = 0.85 
CS7-6: 

Robust Fit : log{L'cs7-&) = 1-03 x \og{Mvir{RcS2-i)) - 1-66 
Least Squares : log{L'cs7-&) = 1.64(±0.05) x \og{Mvir{RcS7-%)) - 3.37(±0.15); r = 0.83 

The robust estimation fit consistently gives values near unity, indicating linear relations be- 
tween L'^^i and Myir- The least squares method yields slight superlinear relations, with values 
around 1.25, except for CS 7-6. More CS 7-6 maps on the lower luminosity may help to clarify 
which fit is better. More importantly, these correlations are by far the strongest we have plotted 
so far, confirming that L'^^^ is an excellent tracer of Myir for the dense gas. 
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In Table [m we give the statistics on the ratios of logarithms of the virial mass to line luminosity 
for each tracer. The scatter is about a factor of 3. All these lines provide good tracers for the mass 
of dense gas, as required for the Gao &: Solomon interpretation. That statement is true even though 
these tracers are optically thick. The explanation is probably similar to explanations given for the 
fact that CO traces the mass of lower density material. Because CO is thermalized at densities 
above about 1000 cm~'^ it is not sensitive to the dense gas that we are studying with these tracers. 
Questions have been raised about the ability of HCN and CS to trace the mass because they require 
high densities to excite (Krumholz & Thompson 2007, Narayanan et al. 2008). In our view, that 
is exactly their virtue: they trace only the dense gas, which CO can fail to trace. The density of 
the material they are tracing depends on their sensitivity to density and temperature, as discussed 
above. What is clear is that they trace the regions where clusters of stars, including massive stars, 
are forming and are thus most relevant to estimating the mass of gas that is immediately relevant 
for star formation. 

The good correlation between L'^^^ and Mvir stands in sharp contrast to the lack of correlation 
between surface density (a good tracer for the mean density of the clump) and L'^^^ (Figures [5T] to 

The spread in L'^^j^ for a given surface density is about three orders of magnitude. We see no 
evidence in the data for the relation between L'^^^ and mean density posited by Narayanan et al. 
(2008). 



5.2. Relations between Lir and L'^^i 

The next question is whether the relations between L/^j and L'^^^ seen in other galaxies (Gao 
& Solomon 2004a) extend down to the scale of individual clumps in our Galaxy. The initial results 
of our survey of HCN 1-0 indicated that they did, except that there is a threshold, below which 
there is not enough dense gas to sample the IMF. Because Ljji is dominated by massive stars, it 
can trace the star formation rate in the same way as galaxies, only with a reasonable sampling of 
the IMF. Even then, it will be subject to large statistical fluctuations from clump to clump, and 
it may underestimate the star formation rate because it takes some time for massive stars to form 
(Krumholz &: Thompson 2007). Therefore, we must expect substantial scatter in this relationship. 

With these caveats in mind, we can ask two questions. Does Ljji correlate with L'^^jl Does the 
relation steepen for tracers with higher critical density, as predicted by Narayanan et al. (2008)? 
We now consider these questions for the four tracers in this survey, as well as CS 5-4 mapping 
survey from Shirley et al. (2003). In Figure E] we plot Lir versus L'^^i for HCN 1-0, HCN 3-2, CS 
2-1, CS 5-4 and CS 7-6. Figure H5] shows the correlation between the distance independent ratio 
Ljr/L'^^i^, an indicator of star formation efficiency, and Ljr for different tracers. Extragalactic data 
are lacking for the CS transitions, so we show the extragalactic data only for the HCN transitions. 
It is important to note that Ljr is the same for each tracer, as it is a property of the entire 
clump and is measured from emission in large beams (e.g., IRAS) or by integrating over maps at 
submillimeter wavelengths. This point has been missed by some analysts. 
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First we consider the plots using L^^^(l — 0) as the dense gas tracer.. The correlation between 
Liji and L^^^(l — 0) found in the galaxies extends to the scale of dense clumps in our Galaxy, as 
originally found by Wu et al. (2005). There is a threshold L'^^^, or equivalently, Lm, which is most 
easily seen in Figure HSl For Ljji > 10^'^ Lq, all dense gas tracers show a nearly linear correlation 
with L/ij. Above the threshold, a linear least squares fit with uncertainties in both axes gives 

logiLji^) = 1.07(±0.06) x log{L'HCNi-o) + 2.98(±0.14); r = 0.85 
and a robust fit gives 

(Lir) = 1.03 X \og{L'Hcm-o) + 3-12. 

The fit to this relation in other galaxies is shown on Figure |3H it is remarkably similar to that for 
massive, dense clumps in our Galaxy. 

Data for other lines in galaxies are scarce. The other panels on Figure 03] and HS] show the 
relations for the clumps in this study, along with the fit for the HON 1-0 line in other galaxies. 
All show correlations with Ljji above the threshold. We give here (in order of increasing ricrit) the 
linear least squares fit with uncertainties in both axes, and the robust fit for other dense gas tracers 
for cores with L//j > 10^'^ Lq. 

CS 2-1: 

Least Squares : log{LiR) = 1.03(±0.05) x log(L^52-i) + 3-25(±0.11); r = 0.80 
Robust Fit : log{LjR) = 0.87 x log{L'cs2~i) + 3.56 

CS 5-4: 

Least Squares fit : log{LiR) = 1.05(±0.05) x log(L^55_4) + 3.77(±0.08); r = 0.86 
Robust Fit : log{Lin) = 0.86 x log{L'cs5-4) + 3.90 

CS 7-6: 

Least Squares fit : logiLm) = 0.81(±0.04) x log(L^57_6) + 4.31(±0.06); r = 0.81 
Robust Fit : log{LiR) = 0.64 x log(L^s7_6) + 4.58 

HCN 3-2: 

Least Squares fit : log{LjR) = 0.88(±0.06) x log{L'jjcN3-2) + 3.94(±0.11); r = 0.82 
Robust Fit : log{LjR) = 0.79 x log(L'^c'^3_2) + 4.09 
All the relations are close to linear. There is some tendency for the tracers with higher Ticrit 

to 

have slopes less than unity, as predicted by Narayanan et al. (2008). However they predict a slope 



- 25 - 



of about 0.7 (ranges from 0.63~0.73 from their Fig. 7) for HCN 3-2 whereas we see 0.88 it 0.06. 
They do not predict values for CS transitions. 

The ratio of Ljji to L'^^^ depends on the tracer. The constant term in the fit increases steadily 
as a function of excitation requirements. As discussed in §4.1, the luminosity of the higher J tracers 
is less for a given clump because of their higher requirements for both density and temperature. It 
is then natural that the relation depends on the tracer used. The dependence on critical density 
(or more generally on excitation requirements) is primarily in this constant term (which increases 
by an order of magnitude from CS 2-1 to CS 7-6) rather than in the slope in the log-log plot. 



5.3. Comparing Explanations 

In order to explain the threshold and the linear correlation over a large range of scales, we have 
proposed a model of clustered star formation (Wu et al. 2005). We suggested that there is a basic 
unit of clustered formation. For Mdense less than the mass of this unit, Ljji/Mdense rises rapidly 
with Mdense, 8tS higher mass stars can form, and stellar luminosity rises rapidly with stellar mass. 
For Mdense greater than the mass of this unit, the IMF is reasonably sampled and further increases 
in mass produce more units, but no further change in Ljji/Mdense- If we suppose that larger 
scale cluster formation is built up by adding more and more such units, then the linear correlation 
between the total Lj^ and Mdense is a natural result, since L'^^j^ traces Mdense- In that case, 
the only difference between star formation on different scales and in different environments- from 
molecular clouds in our Galaxy to massive starbursts- is just how many such units they contain. 
Table [TSl and Table [TBI contain the statistics of infrared luminosity, virial mass, and L/ij/M^j^mi of 
all the massive dense clumps, and of the sample of massive dense clumps with Lm > 10^'^ Lq only, 
respectively. 

The advantage of this model is that it can easily explain not only the linear correlation, but 
also the luminosity threshold: above a luminosity of about 10^ Lq, the IMF is reasonably well 
sampled, and the dependence of stellar luminosity on stellar mass is less steep. The luminosity 
threshold is set by the requirement to sample the IMF. A threshold of > 10^'^ Lq is consistent 
with a value (~ IO^Lq) that dominates the luminosity function of the UCHII regions in the Milky 
Way (Cassassus et al. 2000). This model can explain the fact that starburst galaxies have a very 
high star formation rate, but a similar star formation efficiency with respect to the dense gas as 
individual massive dense clumps in the Milky Way. Extreme starbursts in galaxies happen when 
entire molecular clouds become filled with dense gas, likely due to the huge pressure in the extreme 
environments. 

To summarize, the model proposed by Wu et al. (2005) starts with what we know about 
well-studied regions of star formation in our Galaxy and tries to extend that picture to other 
galaxies. 

The alternative models assume that star formation laws derived from extragalactic studies cov- 
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ering large scales can be extended to the scales where stars actually form. For example, Krumholz 
& Thompson (2007) assumed that a Schmidt law in the form involving volume density describes 
star formation: oc p]^'^ over the full range of densities. They further assumed that the distri- 
bution of densities in a molecular cloud follows a lognormal probability distribution. They argued 
that a molecular tracer essentially measures the mass of gas with density equal to or greater than 
the critical density for that transition. With a simple model of radiative transfer to correct the 
critical density for trapping, they showed that these assumptions led to a linear relation between 
star formation rate and line luminosity of a transition when the mean density of the galaxy was 
less than the critical density of the transition, as seen for HCN 1-0. If, on the other hand, the mean 
density of the gas ((n)) is greater than the critical density of the tracer, then the Lm — L[^acer 
correlation will be super-linear, as is the case for CO. Narayanan et al. (2008) add to this picture 
the idea that index is flattened by the increasing index (/3) in the relation between L'^^^ and density. 

The average density determined from CS excitation of the massive clumps in our sample is 
about 10^'^ cm~'^ (Plume et al. 1997), less than the critical density of all the tracers in this study 
except for the CS 2-1 line (Table [9]), but greater than the effective density (Table [9]) and the 
density that was found to contribute most to the HCN 1-0 line in the simulations of Krumholz and 
Thompson (2007). In fact, a density derived from excitation analysis is biased towards the densest 
regions and the mean density of the clumps in the sense of mass divided by volume is generally less 
(e.g., Shirley et al. 2003). As noted above, the relations we find do not support the suggestions by 
Krumholz & Thompson (2007) or Narayanan et al. (2008). 

Another problem with the Krumholz-Thompson argument is that there is no evidence in our 
Galaxy that the star formation density depends on the local gas density in the way assumed 
for a Schmidt law. Rather, there is a threshold density below which there is essentially no star 
formation. Above that density, clumps and cores do not have lognormal distributions, but instead 
have centrally condensed power law distributions (Mueller et al. 2002, Shirley et al. 2003). 

More detailed discussion about the origin of the Ljr — L[,^acer correlation is beyond the scope 
of this paper. More discussion can be found in recent reviews and papers (e.g., Shirley et al. 2008; 
Narayanan et al. 2008; Evans 2008; Baan et al. 2008; Inno et al. 2009; Juneau et al. 2009). 

5.4. Caveats: Different Beams and Different Distances 

The beam sizes used to map the lower excitation lines were roughly twice as large as those used 
for the higher excitation lines (Table 2). As shown by Liszt & Linke (1975), derivation of densities 
via excitation analysis with such data can be very misleading. Densities from such an analysis with 
better-matched beams have been obtained for this sample by Plume et al. (1997), and we have not 
analyzed the current data in that way. Instead, we focused on the new feature of these data, which 
are maps of multiple transitions. We treated each transition separately and derived mean surface 
and volume densities for each tracer (§4.4). We then discussed the difference between these mean 
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densities and those derived by Plume et al. (1997) from excitation analysis. 

The sources in this study lie at a large range of distances (0.2 to 13.7 kpc). This is an inevitable 
consequence of our attempt to capture the full range of properties of massive dense clumps, as many 
of the most extreme star- forming clumps are at large distances. The effective linear resolution is 
much worse for the distant clumps, which could bias the determinations of some quantities, such as 
surface and volume densities, to lower values. Examination of equations 3,5, and 6 would suggest 
that S cx and n oc d~'^ if there were no trends in mass with distance. To check the importance of 
this effect, we plot in Figure [M] and Figure HT] the logarithms of surface and mean volume densities 
versus the logarithm of distance. The lines in the plots have slopes fixed to the values suggested 
above and only the offset is fitted. With substantial scatter, the expected relations are seen. Note 
that these relations assume that there is no distance effect on the clump masses. Sensitivity limits 
might lead to more distant clouds being systematically more massive, partly canceling the effects 
of decreasing resolution. 

We would expect no strong trends with distance in line ratio plots, and indeed none are seen 
(Figure 08]). 

6. Summary 

Maps of over 50 massive, dense clumps in HCN and CS were obtained. The main sample was 
drawn from the sample of Plume et al. (1992) selected from water maser surveys and found to have 
CS 7-6 emission. In addition, we mapped a smaller number of sources of lower luminosity from other 
surveys. These are not included in statistics quoted for massive, dense clumps. We have presented 
maps of HCN 1-0 and 3-2, along with CS 2-1 and 7-6. Together with previously published maps of 
CS 5-4, these provide the basis for analysis of systematics of massive, dense clumps in tracers with 
a wide range of excitation requirements. Upper state energies in temperature units range from 4.3 
K to 66 K, critical densities range from 3.9 x 10^ cm~'^ to 6.8 x 10^ cm~'^, and effective densities 
range from 4.1 x 10^ cm~^ to 2.6 x 10^ cm^'^. 

As excitation requirements increase, FWHM map sizes shrink, as expected for centrally peaked 
density distributions. Interestingly, the linewidths increase with excitation requirement, producing 
an inverse size-linewidth relation, opposite to the one found for low density regions. The linewidths 
are highly supersonic, with Mach numbers greater than 6. These facts suggest that outfiows and 
expanding HII regions are injecting turbulence on small scales. An alternative would be increased 
inflow velocities. 

Virial masses were computed for each of the tracers. Because the tracers with higher excitation 
requirements have smaller FWHM sizes, the virial masses derived from them are less. For the tracers 
that best trace the whole clump, CS 2-1 and HCN 1-0, the mean and median virial masses are 5000 
and 1950 Mq (CS), or 5300 and 2700 Mq (HCN). The mean and median infrared luminosities are 
4.7 X 10^ Lq and 1.06 x 10^ Lq. These clumps include some of the most massive and luminous star 
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forming regions in the Galaxy, such as W43, W51, and W49. For clumps with luminosity above 
lO^-^ L0, {Lir) = (6.3 ± 1.6) X 10^ L0. 

Mass surface densities are similar to those needed for high-mass star formation, according to 

theoretical arguments (Krumholz & McKee 2008). Even though only CS 7-6 indicates a mean 
surface density exceeding their threshold of 1 g cm~^ for massive star formation, the observational 
and theoretical uncertainties are sufficient to suggest that most of these clumps can form massive 
stars. And their high luminosities indicate that indeed they are doing so. 

Line profiles skewed to the blue, using the blue peak criterion, are common (16 examples) 
in the sample, but some (9) are skewed to the red. This excess of blue over red is less than in 
the original sample of Wu k, Evans (2003), but it still suggests the existence of large scale inflow 
motions in many clumps. 

The line luminosities all correlate very well with virial mass, with different offsets because of 
the fact that higher excitation transitions trace smaller fractions of the core. There is no apparent 
correlation of line luminosity with surface density or volume density, contrary to some recent sug- 
gestions (e.g., Narayanan et al. 2008). Line ratios show little correlation with any other properties. 
In principle, ratios of line luminosity should trace density, but in a power-law density distribution, 
densities range widely and line ratios tend to constant values. 

The far-infrared luminosity correlates very well with all the line luminosities. The correlations 
are close to linear, with some marginal evidence for lower exponents for lines with the most stringent 
excitation requirements. While that agrees qualitatively with predictions by Narayanan et al. 
(2008), the effect is much less than the predicted one. 

Overall, the data support the idea that the linear relation between Ljr and L'^^^ for starburst 
galaxies and molecular lines that are hard to excite can be explained if such galaxies have large 
numbers of massive, dense clumps similar to those in this study. In this case, the star formation 

rate in starburst galaxies depends on the total mass of gas above the threshold density or surface 
density needed for massive star formation, rather than on a local, non-linear form of the Schmidt 
law. 

We thank D. Narayanan and M. Krumholz for lively discussions that illuminated their mod- 
els. Jingwen Wu thanks the postdoctoral fellowship support from the Submillimeter Array of the 
Smithsonian Observatory while completing this work. This work was supported in part by NSF 
Grant AST-0607793 and by McDonald Observatory. 
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Table 1. Mapping Surveys towards Massive Clumps 



Source a (2000.0) 5 (2000.0) dist ref. ^ CS CS"^ CS HCN HON 350 /im ^ 





















kpc 
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5—^4 


7 — > 6 


1 — > 


3 — > 2 


bHAHG 


G121.30+0.66 


00 


36 
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.5 


63 


29 


02, 


,0 


1.2 


1 


M 


M 


M 


M 


M 
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G123.07-6.31 


00 


52 


25, 


,2 


56 


33 


53, 


,2 


2.2 


2 


M 


M 


M 


M 


M 


M 




02 


25 
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.6 


62 


05 


51, 


,1 


2.4 


18 












M 


W3(OH) 


02 


27 
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.7 


61 


52 


25. 


.5 
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M 
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27 
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.5 


58 


46 
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M 
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05 


39 


12, 


.9 


35 


45 
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.0 
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M 


M 


M 


M 
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M 
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35 
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20 
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Tablo 1 Contimiod 



Source 




(^zOUO.UJ 
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ret. 


Ob 


















kpc 
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7 — > 6 


1 — > 


3—^2 


bHARC 


G48.61+0.02 


19 20 31.4 


13 55 24.2 


11.8 


1 


M 


M 


M 


M 


M 


M 


W51W 


19 23 11.0 


14 26 38.1 


7.0 


17 




M 


M 








W51M 


19 23 43.9 


14 30 29.4 


7.0 


17 


M 


M 


M 


M 


M 




G59.78+0.06 


19 43 11.6 


23 43 53.9 


2.2 


1 


M 


M 




M 






S87 


19 46 20.5 


24 35 34.4 


1.9 


15 


M 






M 




M 


S88B 


19 46 47.7 


25 12 56.2 


2.1 


2 


M 


M 


M 


M 


M 


M 


K3-50 


20 01 45.7 


33 32 41.1 


9.0 


4 


M 


M 




M 






ON 1 


20 10 09.1 


31 31 37.3 


6.0 


1 


M 


M 


M 


M 




M 


ON 2S 


20 21 41.0 


37 25 29.4 


5.5 


1 


M 


M 


M 


M 


M 


M 


ON 2N 


20 21 43.9 


37 26 38.6 


5.5 


1 




M 


M 








S106 


20 27 25.7 


37 22 51.8 


4.1 


16 


M 


M 


M 


M 


M 


M 


W75N 


20 38 36.9 


42 37 37.5 


3.0 


4 


M 


M 


M 


M 


M 




DR21 S 


20 39 00.8 


42 19 29.8 


3.0 


4 


M 


M 


M 


M 


M 




W75(OH) 


20 39 01.0 


42 22 49.8 


3.0 


4 


M 


M 


M 


M 


M 




G97.53+3.19 


21 32 11.5 


55 53 54.6 


7.9 


1 




M 






M 


M 


BFS 11-B 


21 43 06.7 


66 07 04.0 


2.0 


5 




M 


M 






M 


CepA 


22 56 18.1 


62 01 46.2 


0.73 


17 


M 


M 




M 


M 


M 


S158 


23 13 44.9 


61 26 50.6 


2.8 


2 


M 


M 




M 


M 


M 


S158A 


23 13 45.4 


61 28 11.7 


2.8 


2 


M 


M 




M 


M 


M 


NGC7538-IRS9 


23 14 01.7 


61 27 19.9 


2.8 


2 


M 






M 


M 


M 


S157 


23 16 04.4 


60 01 40.6 


2.5 


2 


M 


M 


M 


M 


M 


M 


IRAS19094+0944 


19 11 52.0 


09 49 46.2 


6.11 


19 


M 






M 








90 07 Ofi 7 


97 9K ^9 Q 


0.73 


19 


M 






M 


M 




IRAS20106+3545 


20 12 31.3 


35 54 46.3 


1.64 


19 


M 






M 






IRAS20126+4104 


20 14 25.1 


41 13 32.4 


4.18 


19 


M 






M 






IRAS20188+3928 


20 20 39.4 


39 37 51.9 


0.31 


19 


M 






M 


M 




IRAS20216+4107 


20 23 24.2 


41 17 39.8 


1.70 


20 


M 






M 


M 




IRAS20220+3728 


20 23 55.7 


37 38 10.5 


4.49 


19 


M 






M 






IRAS20333+4102 


20 35 09.6 


41 13 17.7 


1.6 


19 


M 






M 


M 




IRAS21391+5802 


21 40 42.4 


58 16 09.6 


0.75 


19 


M 






M 


M 




IRAS22172+5549 


22 19 09.0 


56 04 44.8 


2.87 


19 


M 






M 






IRAS22198+6336 


22 21 27.6 


63 51 42.3 


1.28 


19 


M 






M 






IRAS2301H-6126 


23 03 13.0 


61 42 26.3 


0.73 


20 


M 






M 


M 




IRAS23385+6053 


23 40 53.2 


61 10 21.0 


6.94 


19 


M 






M 






IRAS03282+3035 


03 31 21.0 


30 45 32.6 


0.22 


21 


M 






M 


M 




L483 


18 17 29.8 


-04 39 38.3 


0.20 


22 










M 




L1251B 


22 38 47.1 


75 11 28.8 


0.30 


23 


M 






M 







* "M" in the table indicates that the source has been mapped in this transition 

t-l. Plume ct al. 1992, 2. Blitz ct al. 1982, 3. Harris & Wynn- Williams 1976, 4.Gcnzcl & Dowries 1977, 5. 
Braz & Epchtein 1983, 6. Chini et al. 1986, 7. Hofner et al. 1996, 8. Val'tts et al. 2000, 9. Solomon et al. 1987, 
10. Zinchecnko et al. 1994, 11. Hunter et al. 2000, 12. Churchwell et al. 1990, 13. Downes et al. 1980, 14. Wood 
& Churchwell 1989, 15. Brand & Blitz 1993, 16. Wink et al. 1982, 17. Zhou et al. 1996, 18. Plume et al. 1997, 19. 
Zhang et al. 2005, 20. Wu et al. 2004, 21. Cernis 1990, 22. Dame & Thaddeus 1985, 23. Kun & Prusti 1993. 

'^Mapping survey done by Shirley et al. (2003). 



■^Mapping survey done by Mueller et al. (2002). 
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Table 2: Observing Information 



iransitions 


Dates 


ieiescope 


rrequency 
(GHz) 


/) 

"rah 

(arcsccj 


rimb 


HCN 3-2 


Jun 2002 


CSO 


265886.188 


28.1 


0.64 




Dec 2002 


CSO 


r\ r- o o -1 act 

265886.188 


28.1 


0.60 




Jul 2003 


CSO 


265886.434 


28.1 


0.76 




Jan 2004 


CSO 


265886.434 


28.1 


0.66 




Jun 2004 


CSO 


265886.434 


28.1 


0.64 




Oct 2007 


CSO 


265886.434 


28.1 


0.77 


CS 5-4 


Sep 1996 


CSO 


244935.610 


24.5 


0.56 




Dec 1996 


CSO 


244935.610 


24.5 


0.54 




Apr 1997 


CSO 


244935.610 


24.5 


0.56 




Jun 1997 


CSO 


244935.610 


24.5 


0.58 




Dec 1997 


CSO 


244935.610 


24.5 


0.56 




Jul 1998 


CSO 


244935.610 


24.5 


0.56 




Jul 1999 


CSO 


244935.610 


24.5 


0.60 


CS 7-6 


Jul 1998 


CSO 


342882.900 


16.9 


0.59 




Dec 1998 


CSO 


342882.900 


20.5 


0.61 




Jun 2000 


CSO 


342882.900 


24.7 


0.77 




Dec 2000 


CSO 


342882.900 


24.7 


0.63 




Jun 2002 


CSO 


342882.900 


24.7 


0.53 


HCN 1-0 


Apr 2004 


FCRAO 


88631.847 


58.0 


0.50 




Dec 2004 


FCRAO 


88631.847 


58.0 


0.50 




Feb 2005 


FCRAO 


88631.847 


58.0 


0.50 


CS 2-1 


Apr 2004 


FCRAO 


97980.968 


50.0 


0.50 




JJec zuU4 


T^r^R A Pi 
r UxlAU 


y ( you.yoo 


oU.U 


u.ou 




Feb 2005 


FCRAO 


97980.968 


50.0 


0.50 


C^^S 5-4 


Dec 1998 


CSO 


241016.176 


31.0 


0.66 




Jul 2001 


CSO 


241016.176 


31.0 


0.73 




Jan 2002 


CSO 


241016.176 


31.0 


0.53 




Jun 2002 


CSO 


241016.176 


31.0 


0.59 




Nov 2005 


CSO 


241016.089 


31.0 


0.80 


H^^CN 3-2 


Dec 2002 


CSO 


259011.814 


28.9 


0.60 




May 2003 


CSO 


259011.814 


28.9 


0.58 




Jun 2004 


CSO 


259011.814 


28.9 


0.64 




Jun 2006 


CSO 


259011.814 


28.9 


0.80 



Table 3. Observational Results of CS and C^^ Transitions 



SOURCE 






CS 2 ^ 1 






CS 7 


6 






C^^S 5 -)■ 4 






T* 

R 




^LSR 


FWHM 


R 




^ LSR 


FWHM 


T* 

R 






FWHM 




(K) 


(Kkms 




(kms-1) 


(kms~^) 


(K) 


(Kkms~^) 


(kms~^) 


(kms~^) 


(K) 


(Kkms-l) 


(kms-i) 


(kms~^) 


G121.30+0.66 


4.2(0.2) 


18.9 


0.8) 


-17.87(03) 


3.82(09) 


3.1(0.4) 


9.0(1.1) 


-17.44(05) 


2.58(13) 


0.2(.08) 


0.6(0.2) 


-17.64(17) 


2.87(37) 


G123.07-6.31 


3.6(0.3) 


13.9 


1.0) 


-30.56(04) 


3.56(10) 


2.6(0.5) 


14.3(1.6) 


-30.97(09) 


5.33(21) 


0.2(.15) 


1.3(0.4) 


-30.29(37) 


5.35(108) 


G135.28+2.80 












0.9(0.3) 


2.8(0.8) 


-71.49(09) 


2.47(26) 










W3(OH) 


6.7(0.1) 


36.0 


0.6) 


-46.80(01) 


5.02(03) 


8.1(0.7) 


50.9(2.4) 


-47.65(04) 


6.05(10) 


1.0(0.1) 


5.9(0.5) 


-47.49(08) 


5.50(20) 


S231 


3.6(0.4) 


18.4 


1.5) 


-16.74(06) 


4.76(16) 


3.5(0.4) 


17.3(1.5) 


-15.82(04) 


4.38(11) 


0.5(.08) 


1.2(0.2) 


-15.77(07) 


2.48(15) 


S235 


4.1(0.2) 


10.9 


0.8) 


-16.78(02) 


2.33(06) 


5.0(0.2) 


19.5(0.8) 


-17.09(02) 


3.38(04) 


0.8(0.2) 


2.0(0.5) 


-16.93(09) 


2.35(24) 


S241 


2.5(0.1) 


7.3 


'0.4) 


-8.74(02) 


2.63(05) 


0.6(0.2) 


1.9(0.5) 


-8.79(11) 


2.98(24) 










S252A 


3.8(0.2) 


12.8 


0.8) 


9.05(03) 


3.09(07) 


2.7(0.3) 


10.5(0.9) 


9.63(05) 


3.57(11) 










S255 












7.2(0.4) 


32.6(1.2) 


7.52(02) 


4.03(05) 


1.1(0.2) 


3.9(0.7) 


7.34(09) 


3.03(25) 


RCW142 


7.3(0.3) 


44.2 


1.4) 


17.08(03) 


5.23(08) 


10.6(0.5) 


83.0(2.3) 


16.63(03) 


6.81(07) 


4.3(0.4) 


26.5(1.4) 


17.35(05) 


5.60(12) 


W28A2(1) 


5.4(0.3) 


30.1 


1.1) 


9.38(03) 


4.77(09) 


12.5(0.4) 


131.5(2.4) 


10.40(03) 


9.29(09 


3.4(0.2) 


23.4(0.7) 


9.77(03) 


5.91(09) 


MS E 


4.3(0.3) 


10.4 


'0.8) 


10.58(03) 


2.15(08) 


5.3(0.4) 


18.4(1.2) 


10.91(03) 


3.29(07) 


1.9(0.2) 


4.6(0.4) 


10.81(03) 


2.23(08) 


G9.62+0.10 


2.3(0.5) 


17.9 


2.4) 


4.72(17) 


7.27(43) 


5.9(0.4) 


48.9(2.2) 


5.42(04) 


7.42(11) 


1.8(0.1) 


15.7(0.8) 


4.28(06) 


6.73(17) 


G8.67-0.36 


1.7(0.3) 


14.3 


1.1) 


35.23(13) 


8.12(32) 


4.6(0.3) 


30.1(1.5) 


34.97(04) 


5.74(10) 


1.0(0.2) 


5.7(0.5) 


34.84(08) 


5.08(23fiL 


W31 (1) 












3.0(0.9) 


19.9(3.2) 


65.00(18) 


5.84(43) 


0.9(0.1) 


8.0(0.6) 


66.27(12) 


8.56(29r 


GlO.60-0.40 


4.2(0.3) 


31.3 


1.5) 


-3.24(06) 


6.67(14) 


13.0(0.5) 


131.3(2.5) 


-3.56(03) 


9.17(07) 


4.1(0.2) 


29.8(0.8) 


-2.86(03) 


6.72(08) 


G12.42+0.50 


2.4(0.3) 


11.4 


1.1) 


18.18(07) 


4.32(18) 


















G12.89+0.49 


1.5(0.4) 


8.1 


1.6) 


33.77(16) 


4.39(44) 


3.9(0.3) 


21.1(1.1) 


33.18(04) 


5.09(09) 


1.3(0.1) 


5.2(0.3) 


33.22(03) 


3.78(08) 


G12.21-0.10 


1.0(0.2) 


7.3 


'0.7) 


24.92(12) 


7.16(28) 


2.4(0.3) 


18.8(1.5) 


23.72(08) 


7.63(19) 


0.6(.07) 


4.1(0.3) 


24.44(08) 


7.06(18) 


W33cont 












8.4(0.4) 


69.5(1.5) 


35.45(02) 


7.79(06) 


3.8(0.1) 


21.2(0.5) 


35.25(02) 


5.13(05) 


G13.87+0.28 


11.2(0.4) 


76.3 


1.8) 




6.38(06) 


4.6(0.4) 


14.4(1.1) 


48.32(03) 


2.57(07) 


0.9(0.1) 


2.5(0.4) 


48.57(07) 


2.50(17) 


W33A 


1.4(0.3) 


7.8 


1.4) 


36.58(16) 


5.59(36) 


1.6(0.5) 


10.3(1.4) 


37.59(14) 


5.78(41) 


0.7(0.1) 


2.4(0.3) 


37.89(08) 


3.22(24) 


G14.33-0.64 


1.7(0.3) 


10.7 


1.2) 


22.66(12) 


5.51(33) 


6.6(0.6) 


34.7(2.3) 


22.16(04) 


4.78(11) 


1.5(0.1) 


4.4(0.3) 


22.19(03) 


2.74(08) 


G19.61-0.23 


0.6(0.2) 


4.0 


'0.7) 


42.67(24) 


5.61(62) 


3.1(0.8) 


31.1(3.5) 


42.63(17) 


9.72(39) 










G20.08-0.13 












3.0(0.3) 


29.1(1.6) 


41.74(06) 


8.83(17) 










G23.95+0.16 


3.4(0.2) 


10.4 


'0.8) 


79.54(03) 


2.65(09) 


2.9(0.3) 


8.5(0.9) 


80.57(05) 


2.74(11) 


1.2(0.3) 


2.9(0.6) 


80.26(09) 


2.39(22) 


G24.49-0.04 


1.3(0.2) 


7.3 


0.7) 


110.17(08) 


4.57(23) 


1.3(0.3) 


4.4(1.0) 


109.60(15) 


6.14(38) 










W42 












4.3(0.5) 


36.5(2.3) 


110.90(07) 


7.65(18) 










W43S 


5.5(0.1) 


28.2(0.7) 


97.51(02) 


4.38(05) 


7.2(0.7) 


41.0(2.6) 


97.76(05) 


5.13(11) 


1.4(0.1) 


5.9(0.1) 


98.00(02) 


3.97(04) 


G31.41+0.31 












4.6(0.3) 


37.4(1.4) 


96.22(04) 


7.27(11) 


1.2(0.1) 


7.8(0.5) 


96.97(06) 


5.86(15) 


W43 Main3 












2.3(0.4) 


18.8(1.8) 


98.36(10) 


7.53(25) 










G31.44-0.26 












2.6(0.2) 


12.9(0.8) 


86.75(04) 


4.68(09) 


0.6(.14) 


2.2(0.4) 


87.07(16) 


3.80(36) 


G32.05+0.06 












1.5(0.3) 


11.7(1.1) 


95.50(10) 


7.26(22) 


0.3(.ll) 


1.2(0.3) 


95.71(22) 


6.10(66) 



Table 3 — Continued 



SOURCE 




CS 2 


1 






CS 7 


6 






C^^S 5^4 






R 






FWHM 








FWHM 


R 




^LSR 


FWHM 




(K) 


(Kkms-i) 


(kms-i) 


(kms-i) 


(K) 


(Kkms-i) 


(kms^^) 


(kms~^) 


(K) 


(Kkms-i) 


(kms-^) 


(kms-i) 


G32.80+0.20 










4.4(0.6) 


38.8(2.7) 


14.22(08) 


8.40(18) 










W44 


7.0(0.2) 


32.8(0.7) 


57.62(02) 


4.12(04) 


11.8(0.4) 


74.9(1.9) 


57.49(02) 


5.68(06) 


4.3(0.2) 


Zz.i[\j.D ) 


DO.Z4i[y)Z ) 


4.03(08)* 


G35. 58-0.03 


0.8(0.2) 


4.5(0.6) 


53.22(12) 


5.11(29) 


3.0(0.5) 


19.7(2.3) 


53.22(08) 


5.72(23) 


0.3(.16) 


Z.\J[U. I ) 


oo.\)o[oi- ) 


5.88(82) 


G35.20-0.74 










3.6(0.8) 


28.1(3.5) 


34.69(13) 


7.12(30) 


0.5(0.1) 


o /I /n q\ 
Z.4i[\j.o ) 


oo. t l[ W ) 


4.70(26) 


W49N 


2.1(0.2) 


40.2(0.9)* 


2.54(11)* 


7.51(23) 


8.9(0.6)* 


130.7(3.2) 


6.81(05) 


13.88(11) 


0.5(0.1) 


o.D(^U.4 ) 


Q 1 Q /I 

o. iol^lo ) 


6.41(39) 




2.8(0.2) 




12.00(08) 


7.78(19) 










0.6(0.1) 


A Tff^ A \ 

4. / \y-^) 


lZ.OO[l.t) ) 


7.86(43) 


W49S 










2.8(0.4) 


23.2(1.5) 


15.42(08) 


7.78(18) 


0.2(.09) 




1 A QO/'OQ\ 

J.4:.OZyZO ) 


7.56(54) 


G45.07+0.13 


2.0(0.2) 


12.2(1.1) 


59.05(07) 


5.13(21) 


5.7(0.6) 


43.1(2.8) 


59.28(06) 


7.00(15) 


1.2(0.2) 


Q '7^^\ (i\ 

O. { [KJ.D ) 


oy.Zi [yjo ) 


6.33(23) 


G48.61+0.02 


2.0(0.1) 


9.3(0.5) 


18.45(04) 


4.08(11) 


1.8(0.3) 


12.0(1.4) 


18.00(08) 


6.00(24) 


0.2(.06) 


n c;/ 1 c;\ 
U.o^. io ) 


1Q 1Q/'^1^ 
io. iol^li ) 


2.34(21) 


W51W 










5.2(0.6) 


23.6(2.0) 


51.04(05) 


4.12(12) 


0.9(0.2) 




Ol.UH^Uo ) 


3.39(25) 


W51M 


6.3(0.2) 


71.8(1.1) 


56.82(03) 


10.57(06) 


13.1(0.4) 


149.6(2.3) 


56.29(02) 


10.69(06) 


2.8(0.3) 


07/1/1 Q 

Z i .4(^1.0 ) 


oo.yy^u i ) 


8.96(17) 


G59.78+0.06 


2.6(0.2) 


7.2(0.7) 


22.78(04) 


2.36(10) 










0.5(.13) 


l.Z[\J.o ) 


09 7/1 /T 
ZZ. i 4l^lU ) 


2.22(31) 


S87 


4.1(0.2) 


13.9(0.6) 


23.21(02) 


3.17(06) 










0.4(0.1) 


1 1 /n o^ 
L.lyxj.z J 


Zo.ooyJty ) 


2.38(2p) 


S88B 


5.0(0.3) 


15.6(0.9) 


21.60(03) 


2.81(07) 


1.7(0.3) 


4.7(0.7) 


21.88(05) 


2.23(12) 


0.5(.06) 


1 Q/n 1 A 
i.ot^U.i ) 


ZL.OO[{JO ) 


2.36(1S) 


K3-50 


1.3(0.2) 


10.6(0.8) 


-24.75(09) 


7.24(25) 










0.1(.06) 


n Q/'n ')\ 


-zo.uy i^o± ) 


7.61(t@) 


ON 1 


2.5(0.1) 


15.2(0.5) 


11.56(03) 


5.44(08) 


3.1(0.6) 


12.7(1.9) 


11.44(08) 


4.18(18) 


0.5(.13) 


z.iyKj.o ) 


1 1 TTl/l 
il. ( Ui^lo ) 


4.21(2'8) 


ON 2S 


4.1(0.3) 


24.2(1.2) 


-1.13(05) 


5.22(12) 


4.5(0.7) 


22.1(2.3) 


-1.58(06) 


4.54(15) 


0.5(.07) 


9 OiTl 9^ 


-l.Zi}\\JO ) 


3.65(14) 


ON 2N 










3.8(0.7) 


21.4(2.8) 


0.28(09) 


5.23(20) 


0.8(0.1) 


O.Z\\J.O ) 


\J.O^\\JO ) 


3.72(11) 


S106 


2.8(0.1) 


10.2(0.4) 


-1.28(02) 


2.92(07) 


1.8(0.4) 


9.4(1.6) 


-2.37(11) 


4.91(24) 










W75N 


6.5(0.2) 


29.0(0.7) 


9.41(02) 


4.14(04) 


11.3(0.5) 


65.4(2.1) 


9.48(02) 


5.02(06) 


1.5(0.1) 






4.22(11) 


DR21 S 


















1.4(0.2) 


8.1(0.6) 




4.94(19) 


W75(OH) 


6.7(0.2) 


37.4(0.7) 


-3.57(02) 


5.31(04) 


15.6(0.4) 


92.3(1.5) 


-3.64(01) 


5.47(03) 


1.3(0.2) 


7.4(0.7) 


-2.73(08) 


5.45(17) 


BPS 11-B 










1.5(0.4) 


5.8(1.4) 


-9.79(12) 


3.73(31) 










CcpA 


5.1(0.3) 


25.5(1.1) 


-11.05(03) 


4.58(08) 










0.8(.15) 


3.0(0.4) 


-10.39(09) 


3.68(19) 


S158 


7.8(0.3) 


45.7(1.3) 


-56.56(03) 


5.25(06) 










1.2(0.1) 


5.2(0.4) 


-55.74(05) 


4.30(13) 


S158A 


















1.1(0.3) 


4.1(0.9) 


-57.21(11) 


3.64(27) 


S157 


2.9(0.2) 


11.5(0.8) 


-43.57(04) 


3.50(10) 


1.9(0.4) 


7.1(1.1) 


-44.70(07) 


3.53(18) 










IRAS18586+010 


0.9(.15) 


3.8(0.5) 


38.42(10) 


3.60(29) 


















IRAS 19094+0944 


0.8(.14) 


2.1(0.5) 


65.47(08) 


2.72(20) 


















IRAS20050+2720 


1.5(0.3) 


5.5(0.8) 


6.48(09) 


3.19(22) 


















IRAS20106+3545 


1.5(0.1) 


3.4(0.3) 


7.97(03) 


2.04(08) 


















IRAS20126+4104 


2.6(0.2) 


9.3(0.5) 


-3.52(03) 


3.09(07) 


















IRAS20188+3928 


3.6(0.2) 


16.3(0.9) 


2.05(03) 


3.70(09) 










0.2(.12) 


0.7(.35) 


1.74(33) 


3.75(92) 



Table 3 — Continued 



SOURCE 




CS 2 


1 






CS 7 ^ 6 




C^-iS 5^4 








^LSR 


FWHM 




/T|jdi. Vlsr FWHM 




/T|jdz. Vlsr FWHM 




(K) 


(Kkms-i) 


(kms""*") 


(kms~^) 


(K) 


(Kkms-i) (kms-i) (kms"!) 


(K) 


(Kkms-l) (kms-i) (kms"!) 


IRAS20216+4107 


1.1(0.2) 


2.7(0.5) 


-1.70(06) 


2.12(17) 










IRAS20220+3728 


2.0(0.2) 


8.4(0.6) 


-2.16(05) 


3.61(13) 










IRAS20333+4102 


3.5(0.1) 


6.9(0.4) 


8.55(01) 


1.84(03) 










IRAS21391+5802 


3.9(0.2) 


13.1(0.5) 


0.71(02) 


2.99(05) 










IRAS22172+5549 


1.4(0.2) 


4.0(0.6) 


-43.36(07) 


2.61(17) 










IRAS22 198+6336 


2.1(0.1) 


5.1(0.4) 


-11.09(02) 


2.05(06) 










IRAS23011+6126 


1.4(0.2) 


3.8(0.5) 


-11.22(06) 


2.07(18) 










IRAS23385+6053 


1.8(0.2) 


7.1(0.7) 


-49.47(07) 


3.90(14) 











Table 4. Observational Results of HCN 3-2 and R^^CN 3-2 Transitions 



SOURCE H"CN 3 -)■ 2 HCN 3^2 











FWHM 




^peak 






Line" 


5v 




(K) 


(Kkms-l) 


(kms-i) 


(kms-i) 


(K) 


(kms-l) 


(Kkms-l) 


(kms-i) 


Profile 




G121.30+0.66 


0.5(.ll) 


1.8(0.3) 


-17.77(10) 


3.06(25) 


4.1( 


.30) 


-18.42(22) 


22.6(1.1) 


-17.00(22) 


db 


-0.37(08) 


G123.07-6.31 


0.5(.ll) 


2.8(0.3) 


-30.79(13) 


4.87(32) 


2.5( 


.23) 


-32.75(22) 


22.5(1.0) 


-28.90(22) 


db 


-0.46(12) 


W3(2) 


0.7(.12) 


3.3(0.4) 


-38.73(10) 


4.48(25) 


5.9( 


.29) 


-36.97(32) 


32.8(1.4) 


-39.77(32) 


dr 


0.39(04) 


W3(OH) 


1.8(.13) 


10.1(0.5) 


-47.16(05) 


5.30(12) 


6.6( 


.55) 


-49.23(22) 


51.3(2.3) 


-46.15(22) 


db 


-0.32(03) 


G135.28+2.80 










0.7( 


.21) 


-42.18(16) 


3.3(0.7) 


-42.18(19) 






GL490 










2.4( 


.55) 


-13.18(22) 


11.5(1.7) 


-13.18(11) 






S231 


0.8(.10) 


3.3(0.4) 


-16.41(06) 


3.74(17) 


3.5( 


.23) 


-14.12(22) 


20.2(1.0) 




rbs 


0.67(07) 


S235 


l.l(.ll) 


3.4(0.3) 


-17.44(05) 


2.71(13) 


5.0( 


.31) 


-17.36(11) 


29.1(1.0) 


-17.36(03) 




-0.18(06) 


S241 










1.1( 


.21) 


-8.58(22) 


5.4(0.7) 


-8.58(11) 






MonR2 










7.8( 


.23) 


10.19(16) 


39.0(0.9) 


10.31(02) 






S252A 


0.6(.09) 


1.7(0.3) 


9.13(07) 


2.61(16) 


3.7( 


.32) 


10.04(22) 


20.5(1.2) 


9.65(06) 


rbs 


0.35(05) 


S255 










7.9( 


.68) 


7.98(22) 


48.1(2.5) 


7.65(05) 




0.21(05) 


RCW142 


2.7{.27) 


21.6(1.0) 


16.16(09) 


7.15(23) 


17.8( 


.45) 


14.36(11) 


164.9(2.6) 


20.77(21) 


db 


-0.53(02) 


W28A2(1) 


5.1{.21) 


47.8(1.3) 


9.23(03) 


6.93(09) 


13.1( 


.35) 


11.18(22) 


177.3(2.3) 


11.40(03) 


ws 


0.24(02) 


MS E 


0.8(.13) 


2.5(0.4) 


10.95(06) 


2.76(14) 


7.1( 


.32) 


9.54(22) 


45.6(1.4) 


11.07(10) 


d 




G9. 62+0. 10 


1.2(.19) 


10.9(0.8) 


4.77(10) 


8.17(26) 


6.2( 


.60) 


7.38(22) 


52.0(2.9) 


2.67(42) 


dr 


0.46(02) 


G8. 67-0. 36 


1.4(.33) 


6.9(1.0) 


34.32(13) 


4.40(30) 


2.7( 


.16) 


31.91(16) 


22.4(0.8) 


36.19(50) 


db 


-0.58(04) 


W31 (1) 


l.l(.ll) 


7.9(2.0) 


64.85(08) 


6.69(19) 


6.9( 


.22) 


62.71(11) 


68.2(1.1) 


68.70(21) 


db 


-0.42(02) 


GlO.60-0.40 


3.6(.21) 


30.5(0.9) 


-4.16(06) 


7.74(14) 


11.8( 


.62) 


-6.50(11) 


132.2(3.4) 


-0.43(10) 


db 


-0.54(02) 


G12.42+0.50 


0.9(.26) 


1.7(0.6) 


17.36(12) 


1.79(29) 


4.2( 


.25) 


16.06(11) 


24.1(0.9) 


17.49(09) 


db 


-0.73(15) 


G12. 89+0.49 


0.8(.19) 


5.4(0.8) 


33.13(18) 


6.42(48) 


2.6( 


.15) 




21.1(0.7) 


33.49(30) 


d 




W33cont 


2.1{.34) 


11.6(1.4) 


35.43(11) 


5.74(24) 


10.0( 


.38) 


39.53(22) 


94.3(2.1) 


33.75(20) 


dr 


0.83(02) 


W33A 


0.7(.23) 


4.2(0.8) 


37.69(22) 


5.40(58) 


1.8( 


.21) 


40.49(22) 


25.8(1.2) 


38.14(11) 


rbs 


0.81(07) 


G14.33-0.64 


1.8(.30) 


8.8(1.2) 


22.25(10) 


4.59(25) 


7.1( 


.25) 


25.47(11) 


69.1(1.5) 


21.65(10) 


dr 


1.20(05) 


G19.61-0.23 


0.4(.ll) 


6.3(0.7) 


42.57(34) 


15.50(85) 


3.1( 


.32) 


40.40(22) 


45.2(1.9) 


40.66(14) 


brs 


-0.14(02) 


G20.08-0.13 










2.1( 


.30) 


38.46(32) 


22.3(1.5) 


39.85(16) 


brs 




G23. 9.5+0.16 


0.4{.09) 


1.3(0.3) 


80.32(10) 


3.24(35) 


1.9( 


.28) 


78.57(22) 


14.5(1.1) 


80.52(22) 


db 


-0.71(09) 


W42 


0.9(.32) 


8.4(1.4) 


110.49(28) 


8.55(58) 


3.9( 


.35) 


114.49(22) 


28.1(1.8) 


108.59(22) 


dr 


0.47(05) 


W43S 


1.3(.16) 


7.1(0.7) 


98.11(06) 


5.25(16) 


6.3( 


.32) 


96.61(22) 


52.0(1.4) 


96.85(03) 


brs 


-0.35(03) 


G31.41+0.31 


1.2(.30) 


7.7(1.2) 


96.23(16) 


5.41(46) 


3.2( 


.28) 


93.25(16) 


21.5(1.5) 


98.67(32) 


db 


-0.63(03) 


G31.44-0.26 










3.1( 


.19) 


87.87(16) 


24.8(0.8) 


87.27(05) 




0.21(06) 


G32.05+0.06 


0.4(.16) 


3.8(0.6) 


95.41(28) 


7.17(91) 


3.1( 


.29) 


92.08(16) 


26.2(1.4) 


97.25(32) 


db 


-0.60(08) 


G32.80+0.20 










3.2( 


.25) 




40.8(1.3) 


13.74(99) 


flat 





Table 4 — Continued 



SOURCE Hi^CN 3 2 HCN 3^2 











FWHM 










Line'"' 


6v 




(K) 


(Kkms-i) 


(kms"-*-) 


(kms"-*^) 


(K) 


(kms-l) 


(Kkms-i) 


(kms-1) 


Profile 






Q Q/ ^ A^ 
0.0^^.14 ) 


Z4.4(_U.0 J 


Qolf\A\ 
OO.c5z(^U4 ) 


K liA /I nA 
0.04(^iU ) 


11.5( 


.32) 


04. / y(^lD ) 


o4.y^^i.c5 ) 


Dl. i / {22 ) 


QD 


n QfifrxA \ 
-KJ.tSoyKJ^) 












7.8( 


.23) 


00.00(^10 ) 


Oo.Dt^i.U ) 


Qo nn/'no^ 
oz.yyi^Uz ) 


rbs 




LxOO.Oo — U.Uo 


U.o(. Ub ) 


1.8(0.2) 


0O.Do(lo ) 


a 1 Q f'Q /I \ 
0.1o(o4) 


2.8( 


.34) 


CO m / oo^ 
OZ.yi(zZ) 


30.0(1.7) 


oz.yi(uyj 




-U.lo^Ub ) 


Q on n 7/1 
Lxoo.^U— U. ( 4 


U. ( \.2X} ) 


D.0(^U.O ) 


Q /I A1 iOA\ 

o4.Ui(z4 J 




7.9( 


.26) 


QV Q/1 / 1 Ci\ 

o ( .o4(iD ) 


t^A Q/I 0\ 

O4.o|^i.z j 


QQ nQ/'QO^ 

oo.yol^oz J 


ar 


n '7'7lr\P.\ 
\j. ( ( \\jK> ) 


W4yiN 


1 O/ QO^ 

i.z^^.oy ) 








8.7( 


.33) 


Q no f'oo\ 
o.(jZ[ZZ ) 


1 /I o /I /o n\ 
14z.4(^z.U j 


iz.44l^zz J 


two 




W4yb 










2.8( 


.16) 


ic; QO/ii^\ 
iD.oZ[lo) 


QT Q /n ^^ 
o7.o(U.yj 


1 c QO /" n^;^ 
Lo.o2{Uo) 




n nT/n/i A 
U.U7(U4j 


t^UA'i on n iQ 
*Jrl4o.cSU— U.io 


n r;/" oi \ 




44. i / (Ji ; 


o.4( ^ It) 


2.2( 


.25) 


Ad cn/Qo^ 
4D.C5U(^OZ ) 


OQ '7/1 o^ 
Zo. / (_±.zj 


AO OK/'QO\ 

4z.yo(^oz ^ 


ox 


n Qi /nc^\ 
U.oH^Uo ) 


Vjr40.U/ -j-U.io 


n 1 fiA 
U.n(^. io ) 


A Q/n '7^ 
4.o(^U. / ] 




0.4VJ(^.l0 } 


4.3( 


.27) 


c^n c^o/oo^ 
D\).OZyZZ ) 


41. £5(^1. ^ 


Kn f^o/n/iA 
oy .Ozi^Uo. J 


rw 


n ciA fr\o\ 
{J.\J^[{J2 ) 


AQ fn in no 




1 AfCi 'i\ 

1.4(^U.o j 


1 Q n/i /"ooA 
io.U4(^zz J 


Q ^ /ft 7^ 

0.00(^0 / ) 


2.3( 


.34) 


1 Q nn/ 1 1 ^ 


1 Q '7/'i /I ^ 

iC5. ( (^i.4 j 


1 Q no/'^o^ 

io.uyi^uy J 




-U.Uzl^U ( ) 


wol w 










2.5( 


.32) 


C^O T Q QO\ 


io.Of^i.z ) 


oi.44(^uy J 


sr 


n Qf^/nfiA 
U.oo(^UD ) 


VV 1 iVl 




o\i.o\Z.o ) 




li.U4^^o( ) 


14.2( 


.44) 


f^'l 7/l/'QO^ 
Oi. / 4l^oZ ) 


T f^Q O/O Q\ 

loo.yi^z.o j 


oo.Doi^yy ) 


brs 


-U.oy i^Uz ) 




n 1 1 \ 


ci A tn o\ 
U.4(^U.z ) 




1 nft/'Qft\ 


2.5( 


.27) 


OI c^o/'oo\ 
zi.ozi^zz ) 


1 n/ 1 n\ 
iu.Ut^i.U j 


oi 70/'nt;A 
21. ( yi^uo ) 


SD 


-(J.OO[W ) 


oo/ 










4.6( 


.28) 


oo '7^^/'oo^ 

ZZ. ( Dl zz ) 


occ n/i n^ 
zO.y^^i.U^ 


zz.D ( [\Jo ) 


brs 


n OA 

-U.z4l^UD ) 


oo!$rS 


n A( C\'7\ 
U.4(^.U/ ; 


1 Q/n OA 


oi QK/'nQ\ 


Q 1 1 /"I Q\ 


2.8( 


.25) 


oo nfi/'ooA 
zz.Uut^zz ) 


1 K Q/n n^ 
io.c$(^U.yj 


oo r\fifr\t^\ 
22.[jO[\JO ) 




n oo/'nc;\ 
{j.2Z[{JO ) 


K3-50 










4.1( 


.26) 


-22.i)l{22) 


40.3(1.3) 


oi TO /n c \ 

-2^.78(05; 


sr 


0.29(05) 


OIV 1 


yJ.oy.Zo ) 


o/'n qA 


1 1 R«( ^ tA 
il.DD^l ( J 


4. ( 4^4Z ) 


3.8( 


.40) 


lo.o i t^ZZ j 


OV O/'l 

Z / .Z|^ i .O j 


1 n 70('oo^ 
iU. i ZyZZ ) 


ar 


n ^ofr\^\ 
\J.oz\\Jo ) 


CIN ZD 


n '7( 1 fi\ 
U. i[.Ld) 




1 t^C\l 1 1 \ 

-i.OUt^ii j 


A Q1 / Qn\ 
4.0i(^oU^ 


8.6( 


.33) 


n nv/^ooA 
-U.U / [22 ) 


ftQ nf' 1 fi\ 

Do.U(^i.D j 


-l.Di)[(JO ) 


rbs 


U.oo(^U4 ) 


1\.T 1\.T 










4.2( 


.40) 


o t; o /oo\ 
i.0o(zz ) 


46. 1(2.0) 


n c; 1 / nnA 

u.oi(uyj 


rbs 


n ftn /n /I A 
0.60(04) 


OlUO 










3.0( 


.35) 


-Z.OO^^ZZ ^ 


on o/'l >i'\ 


-o.\J\J\\J 1 ) 


rbs 




W75N 


1.9(.22) 


11.1(1.0) 


9.67(06) 


5.12(16) 


9.5( 


.27) 


11.74(22) 


92.3(1.7) 


9.46(11) 


dr 


0.50(03) 


DR21 S 


1.7(.19) 


15.2(0.9) 


-2.16(10) 


8.39(26) 


15.8( 


.36) 


-6.02(22) 


217.7(2.3) 


-2.00(11) 


db 


-0.76(04) 


W75(OH) 


1.8(.26) 


10.3(1.0) 


-3.35(13) 


5.26(28) 


11. 1( 


.43) 


-6.27(22) 


80.7(2.2) 


-2.24(22) 


db 


-0.65(03) 


G97.53+3.19 










1.0( 


.32) 


-71.88(22) 


8.3(1.5) 


-71.88(27) 






CepA 










7.9( 


.37) 


-12.33(22) 


71.5(2.0) 


-10.41(22) 


db 


-0.53(05) 


S158 


0.9(.13) 


5.7(0.5) 


-56.29(11) 


6.09(26) 


5.5( 


.32) 


-59.05(22) 


52.0(1.8) 


-55.28(22) 


db 


-0.77(04) 


S158A 










7.0( 


.82) 


-58.43(22) 


37.2(3.7) 


-58.13(10) 


brs 


-0.34(05) 


NGC7538-IRS9 










3.0( 


.33) 


-54.75(22) 


18.7(1.3) 




rbs 




S157 










3.0( 


.82) 


-45.03(22) 


17.8(3.0) 


-44.53(15) 




-0.36(06) 


IRAS20050+2720 










2.3( 


.25) 


5.94(22) 


13.1(0.9) 


5.94(06) 






IRAS20216+1107 










2.0( 


.28) 


-1.31(11) 


7.9(0.9) 


-1.31(06) 






lKAS2:-!01i+()12() 










(),!)( 


.29) 


~l().!)i(22) 


7. 1(1. 0) 


-10.91(21) 







"Code of line profiles: db, double peak with blue strong; dr, double peaks with red strong; d: double peaks with neither one obviously stronger 
flat, the line has a flat top; rbs, blue shoulder with peak skewed to red; brs, red shoulder with peak skewed to blue; ws: wings on both sides; rw 
red wing; sr, peak skews to red; br, peak skews to blue; two, the line has two known velocity components. 
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Tablo 5. Observational Results for HON 1-0 Transition 



Source 


T% (K) 


T* (K) 


T% (K) 

R V / 


fT*j^du (Kkms-^) 


V r *7 R fkms ^ ) 


FWHM (kms-i 


Line" 




(F=l-0) 


fF=2-l) 


(F=l-1)) 


(all hpf Fs) 




(F=2-l) 


profile 


G121.30+0.66 


0.9(0.2) 


4.0(0.2) 




24.9(1.2) 


-18.31(33) 




ab 


G123.07-6.31 


0.9(0.2) 


2.6(0.2) 


1.8(0.2) 


19.6(1.0) 


-31.13(04) 


3.63(05) 




W3(OH) 


1.9(0.1) 






41.8(0.7) 






ab,bl 


S231 


0.7(0.2) 


2.3(0.2) 


1.3(0.2) 


20.1(1.2) 


-17.26(06) 


4.21(07) 


bl 


S235 


1.2(0.2) 


4.3(0.2) 


1.5(0.2) 


19.0(1.0) 


-17.17(02) 


2.37(03) 




S241 


0.3(0.1) 


1.5(0.1) 


0.7(0.1) 


7.1(0.5) 


-8.90(03) 


2.66(04) 




S252A 


1.0(0.2) 


3.5(0.2) 


1.6(0.2) 


21.6(1.1) 


9.05(03) 


3.28(04) 




RCW142 


3.1(0.3) 


5.6(0.3) 


3.1(0.3) 


73.2(1.6) 


16.73(06) 


5.31(07) 


bl 


W28A2(1) 


1.9(0.2) 


4.0(0.2) 


1.0(0.5) 


28.1(1.2) 


8.65(03) 


3.66(04) 




M8E 


0.8(0.3) 


4.0(0.3) 


1.8(0.3) 


17.6(1.3) 


10.46(03) 


2.52(04) 




G9.62-I-0.10 


1.2(0.3) 


1.0(0.3) 


1.1(0.3) 


17.3(1.7) 


5.50(13) 


4.31(20) 




G8.67-0.36 


1.0(0.3) 


1.0(0.3) 




11.4(1.3) 


33.15(33) 






GlO.60-0.40 


0.9(0.6) 


2.3(1.0) 


1.3(0.5) 


46.0(1.7) 


-3.45(01) 




bf,bl 


G12.42-I-0.50 


0.5(0.3) 


1.7(0.3) 


1.3(0.3) 


15.1(1.3) 


17.59(11) 


3.88(11) 


bl 


G12.21-0.10 




0.5(0.2) 




4.7(0.8) 






bl 


G13.87+0.28 


2.9(0.4) 


4.4(0.4) 


2.7(0.4) 


55.3(1.9) 






bl 


G19.61-0.23 




0.2(0.1) 




1.6(0.5) 








G23.95-I-0.16 


0.6(0.2) 


2.4(0.2) 


1.1(0.2) 


12.4(1.0) 


78.90(04) 


2.72(06) 




W43S 


1.6(0.1) 


2.9(0.1) 


1.0(0.1) 


26.4(0.8) 


96.72(03) 


4.03(04) 




W44 


1.7(0.2) 


2.2(0.2) 




12.7(0.8) 


57.07(04) 


2.93(05) 




G35.58-0.03 




0.5(0.1) 




5.2(0.7) 






bl 


W49N 








53.1(1.2) 






bl 


G45.07+0.13 


0.5(0.2) 


0.8(0.2) 


0.2(0.2) 


10.9(1.1) 


59.58(03) 


5.93(23) 


bl 


G48.61+0.02 


0.5(0.1) 


2.1(0.1) 


1.1(0.1) 


17.9(0.7) 


18.36(03) 


4.28(04) 


bl 


W51M 








50.4(1.1) 






bl 


G59.78+0.06 


0.5(0.2) 


2.1(0.2) 


1.0(0.2) 


12.3(0.9) 


22.45(05) 


3.01(06) 




S87 


1.1(0.2) 


4.8(0.2) 


2.1(0.2) 


31.2(0.8) 


22.88(02) 


3.71(02) 


bl2 


S88B 


0.9(0.2) 


3.9(0.2) 


2.3(0.2) 


25.2(0.9) 


21.54(02) 


3.22(03) 


bl2 


K3-50 








24.3(0.8) 






bl 


ON 1 








16.0(0.7) 






bl 


ON 2S 








2.5(0.8) 






ab,bl 


S106 


0.5(0.1) 


2.0(0.4) 


1.0(0.1) 


12.7(0.4) 


-1.58(02) 


3.22(03) 


bl2 


W75N 


0.8(0.1) 


1.4(0.1) 


1.0(0.1) 


19.3(0.5) 


9.75(07) 


5.54(08) 


bl 


DR21 S 


1.2(0.2) 


3.4(0.2) 


1.4(0.2) 


19.4(1.1) 


-3.49(07) 


5.69(09) 


ab,bl 


W75(OH) 








21.7(1.0) 






ab,bl 


CepA 


1.0(0.2) 






30.6(0.9) 






bl,bf 


S158 


3.2(0.3) 


5.3(0.3) 


1.9(0.3) 


57.9(1.5) 


-57.14(04) 


4.76(05) 


bl 


SloT 


0.7(0.2) 


2.7(0.2) 


1.3(0. 2) 


19. 5(0. 9) 


-13.58(01) 


3.79(05) 


bl2 


IRAS18586-I-0106 


0.3(0.2) 


0.7(0.2) 


0.2(0.2) 


4.1(0.7) 


38.34(10) 


3.10(15) 




IRAS19094+0944 


0.2(0.1) 


0.8(0.1) 


0.5(0.1) 


5.1(0.6) 


65.48(07) 


3.13(08) 




IRAS20050+2720 


0.3(0.2) 


0.5(0.2) 


0.3(0.2) 


5.7(1.1) 


6.06(27) 


4.72(25) 


bl2 


IRAS20106+3545 


0.3(0.1) 


1.3(0.1) 


0.6(0.1) 


5.2(0.5) 


7.80(03) 


2.15(05) 




IRAS20126+4104 


0.9(0.1) 


2.9(0.1) 


1.8(0.1) 


19.6(0.6) 


-3.60(02) 


3.14(02) 


bl2 


IRAS20188-I-3928 


0.7(0.2) 


1.9(0.5) 


1.3(0.4) 


19.8(0.8) 


2.49(06) 


4.42(07) 


ab? bl,bf 


IRAS20216+4107 


0.2(0.1) 


1.3(0.1) 


0.5(0.1) 


5.4(0.6) 


-1.57(04) 


2.32(06) 
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Table 5 — Continued 



Source 


T^e (K) 


T|j (K) 


TJj (K) 


/T^di/ (Kkms-i) 


Ylsr (kms 1) 


FWHM (kms-i 


Line" 




(F=l-0) 


(F=2-l) 


(F=l-1)) 


(all hpf Fs) 


(F=2-l) 


(F=2-l) 


profile 


IRAS20220+3728 


0.4(0.2) 


1.8(0.2) 


0.3(0.2) 


8.9(0.8) 


-2.30(05) 


3.09(06) 




IRAS20333+4102 


0.9(0.1) 


4.4(0.1) 


1.8(0.1) 


16.9(0.6) 


8.60(01) 


2.21(02) 




IRAS21391+5802 


0.9(0.1) 


1.8(0.5) 


1.2(0.2) 


15.0(0.6) 


0.50(03) 


3.30(04) 


ab? 


IRAS22172+5549 


0.2(0.1) 


1.2(0.1) 


0.5(0.1) 


6.2(0.7) 


-43.51(06) 


3.05(07) 


bl2 


IRAS22198+6336 


0.5(0.1) 


1.3(0.1) 


0.5(0.1) 


6.0(0.6) 


-11.14(04) 


2.35(05) 




IRAS23011+6126 


0.4(0.1) 


0.9(0.1) 


0.6(0.1) 


6.7(0.6) 


-11.76(07) 


3.25(09) 




IRAS23385+6053 


0.5(0.2) 


1.9(0.2) 


0.8(0.2) 


10.0(0.8) 


-49.95(04) 


2.91(05) 




IRAS03282+3035 


0.8(0.1) 


1.3(0.1) 


0.8(0.1) 


3.3(0.4) 


6.94(02) 


0.92(04) 




L1251B 


0.7(0.1) 


1.5(0.3) 


0.9(0.2) 


5.0(0.3) 


-4.03(02) 


1.56(02) 


ab 



*Code of line profiles: ab, absorption features on the line; bl, three components are blended; bl2: Only F=2-l and F=l-1 are 
blended; bf, the line can not be well fitted by 3 Gaussian components. 
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Table 6. The FWHM Sizes of Maps 



SOURCE 


HON 3 -> 


2 


HON 1 -> 





CS 7 -> 


6 


OS 2 ^- 1 




R(HCN3-2) ^ 


Sdec ^ 


R(HCNl-O) 


Sdec ^ 


R(CS7-6) 


Sdec ^ 




^dec 




(pc) 


(") 


(pc) 


(") 


(pc) 


(") 


VPCJ 


\ > 


G121.30+0.66 


0.09(01) 


41.3 


0.13(01) 


62.4 


0.05(01) 


17.0 


0.16(02) 


76.0 


G123.07-6.31 


0.11(02) 


20.9 


0.41(06) 


76.2 


0.08(01) 


14.2 


0.38(05) 


70.8 


W3(2) 


















W3(OH) 


0.22(02) 


37.1 


0.59(04) 


101.9 


0.12(01) 


20.5 


0.54(04) 


93.3 


G135.28+2.80 


















GL490 


0.05(01) 


24.2 














S231 


0.19(03) 


34.0 


0.43(06) 


77.8 


0.11(02) 


20.1 


0.44(06) 


79.3 


S235 


0.14(02) 


37.3 


0.38(04) 


96.8 


0.10(01) 


24.8 


0.32(03) 


81.7 


S241 


0.22(05) 


19.7 


0.66(08) 


58.2 


0.07(01) 


18.2 


0.59(12) 


52.2 


MonR2 


















S252A 


0.10(01) 


27.5 


0.47(16) 


130.4 






0.45(07) 


122.8 


S255 


0.13(10) 


41.0 






0.10(02) 


32.1 






RCW142 


0.20(02) 


40.5 


0.59(06) 


121.8 


0.11(01) 


23.6 


0.36(03) 


74.3 


W28A2(1) 






0.37(03) 


58.2 


0.10(01) 


15.1 


0.30(02) 


48.4 


M8E 


0.16(01) 


35.8 


0.48(13) 


109.2 


0.08(01) 


18.5 


0.35(06) 


80.1 


G9. 62+0. 10 


0.45(04) 


26.6 


1.13(12) 


66.4 


0.36(06) 


21.3 


1.06(21) 


62.2 


G8. 67-0. 36 


0.38(04) 


35.0 






0.18(03) 


16.4 


0.94(22) 


86.2 


W31 (1) 


0.70(10) 


24.1 






0.82(23) 


28.2 






GlO.60-0.40 


0.38(02) 


24.4 


1.03(07) 


65.4 


0.43(03) 


27.1 


0.93(07) 


59.0 


G12.42+0.50 






0.46(07) 


90.7 






0.51(09) 


99.2 


G12.89+0.49 










0.10(01) 


11.3 * 






G12.21— 0.10 






5.67(177) 


170.6 




U * 


6.59(166) 


198.5 


W33cont 










0.58(04) 


58.7 






G13.87+0.28 






1.75(15) 


163.7 


0.27(05) 


25.5 


1.30(11) 


122.2 


W33A 










0.18(03) 


16.3 






G14.33— 0.64 


0.13(02) 


21.3 






0.16(02) 


25.0 






G19.61— 0.23 


0.21(04) 


21.4 






0.16(02) 


16.4 


1.67(54) 


172.4 


G20.08— 0.13 










0.17(02) 


20.6 






G23.95+0.1D 






1.53(42) 


108.6 


0.37(07) 


26.6 


1.42(45) 


101.2 


G24.49-0.04 










0.17(07) 


19.7 






W42 










0.33(02) 


15.2 






G28.86+0.07 


















W43S 


0.55(06) 


26.7 


1.84(20) 


89.3 


0.25(03) 


12.1 * 


1.42(12) 


69.0 


G31.41+0.31 










0.24(02) 


12.7 






W43 Main3 










0.23(03) 


13.9 






G31.44-0.26 










0.56(09) 


21.6 






G32.05+0.06 










0.53(16) 


25.8 






G32.80+0.20 










0.59(07) 


15.6 






W44 


0.39(03) 


43.0 


1.38(23) 


153.7 


0.25(02) 


28.4 


0.81(05) 


89.8 


S76E 


















G35.58-0.03 






0.63(18) 


74.5 


0.14(03) 


16.7 


0.56(06) 


66.1 


G35.20-0.74 










0.20(03) 


24.4 






W49N 






2.79(18) 


97.4 


1.09(11) 


32.0 


2.10(13) 


73.4 


W49S 










0.62(22) 


18.3 
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Table 6 — Continued 



SOURCE 


HON 3 


• 2 


HON 1 


CS 7 -> 


6 


CS 2 1 




R(HCN3-2) 


^dec 


R(HCNl-O) 


^dec ^ 


R(CS7-6) 




R(CS2-1) ^ 


Sdec ^ 




(pel 


(") 


(pc) 


{") 


(pc) 


(") 


(pc) 


{") 


OH43.80-0.13 


















G45.07+0.13 


0.56(09) 


23.7 


1.34(20) 


57.1 


0.41(07) 


17.5 


0.77(11) 


32.6 


G48.61+0.02 


0.66(13) 


23.2 


3.59(32) 


125.4 


0.25(05) 


8.6 * 


3.09(65) 


107.9 


W51W 










0.53(12) 


31.1 






W51M 


0.80(14) 


47.0 


2.73(49) 


161.0 


0.47(03) 


27.7 


2.32(15) 


136.5 


G59. 78+0.06 






0.46(04) 


87.1 






0.47(06) 


88.4 


S87 






0.48(03) 


103.2 






0.45(04) 


96.6 


S88B 


0.17(03) 


33.6 


0.58(08) 


114.3 


0.14(02) 


26.7 


0.32(03) 


63.2 


K3-50 






1.53(12) 


70.1 






1.44(17) 


66.1 


ON 1 






1.31(11) 


90.0 


0.25(04) 


16.9 


1.15(11) 


79.3 


ON 2S 


0.59(09) 


44.4 


1.09(14) 


81.8 


0.46(05) 


34.8 


1.09(13) 


81.7 


ON 2N 










0.35(06) 


25.9 






S106 


0.43(05) 


42.8 


1.11(10) 


112.0 


0.25(08) 


24.9 


1.27(13) 


127.8 


W75N 


0.27(03) 


36.9 


0.63(04) 


86.4 


0.20(02) 


27.5 


0.56(06) 


77.7 


DR21 S 


0.24(01) 


32.4 


0.75(05) 


102.5 


0.25(02) 


34.2 


0.78(07) 


106.7 


W75(OH) 


0.37(06) 


51.0 


0.87(08) 


120.3 


0.21(02) 


28.8 


0.79(06) 


109.0 


G97.53+3.19 


0.37(04) 


18.0 














BFS 11-B 










0.06(03) 


12.4 






CepA 


0.09(01) 


48.5 


0.23(02) 


130.4 






0.27(03) 


150.3 


S158 


0.46(05) 


68.3 


0.85(06) 


125.9 






0.77(06) 


113.6 


S158A 


0.37(03) 


55.2 














NGC7538-IRS9 


0.42(10) 


62.3 














S157 


0.19(04) 


30.9 


0.45(06) 


74.5 


0.15(04) 


24.2 


0.43(08) 


71.7 


IRAS 19094+0944 






1.48(45) 


100.0 






1.24(66) 


84.0 


IRAS20050+2720 


0.08(02) 


43.1 


0.23(06) 


127.5 






0.22(11) 


123.8 


IRAS20106+3545 






0.24(06) 


60.3 






0.27(07) 


67.1 


IRAS20126+4104 






0.48(04) 


47.3 






0.49(08) 


48.4 


IRAS20188+3928 


0.03(01) 


35.7 


0.07(01) 


98.1 






0.08(01) 


105.5 


IRAS20216+4107 


0.05(02) 


11.6 * 


0.26(03) 


62.4 






0.20(04) 


49.7 


IRAS20220+3728 






1.40(28) 


128.4 






1.41(50) 


129.3 


IRAS20333+4102 


0.28(08) 


71.1 


0.42(03) 


108.6 






0.40(04) 


103.0 


IRAS21391+5802 


0.06(01) 


34.5 


0.14(02) 


77.0 






0.14(01) 


74.3 


IRAS22172+5549 






0.36(06) 


51.2 






0.40(10) 


56.8 


IRAS22198+6336 






0.31(07) 


100.0 






0.33(07) 


106.7 


IRAS2301 1+6126 


0.03(01) 


17.3 


0.17(03) 


98.7 






0.17(04) 


94.7 


IRAS23385+6053 






1.13(17) 


67.3 






1.06(17) 


63.2 


IRAS03282+3035 


0.02(01) 


46.5 


0.10(03) 


187.0 






0.09(03) 


173.8 


L483 


0.02(01) 


32.2 














L1251B 






0.07(01) 


92.1 






0.09(02) 


127.3 


Mean'= 


0.32(20) 


35(12) 


1.13(109) 


100(31) 


0.28(21) 


22(8) 


1.03(112) 


90(31) 


Median 


0.26 


35 


0.71 


97 


0.23 


21 


0.77 


81 



* 



These sources are either unresolved, or derived sizes are less than half of the beam size for this transition. This may cause 
the underestimation of the line luminosity. We will ignore these luminosities in later statistics. 

''Numbers in the parentheses are uncertainties of FWHM size in the unit of 0.01 pc. 
''The angular FWHM size of the core deconvolved from the beam size. 
"^Statistics do not include sources below the horizontal line in this table. 
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Table 7. Luminosities of Massive Clumps 



Source 




L'(HCNl-O) 


L'(HCN3-2) 


L'(CS7-6) 


L'(CS5-4) b 


L'(CS2-1) 




(Lq) 


K kms~^pc^ 


K kms~^pc^ 


K kms~^pc^ 


K kms~^pc^ 


K kms~^pc^ 



G121. 30+0.66 


2.85 


0.53(12) 


0.09(17) 


-0.29(26) 


0.18(13) 


0.45(15) 


G123. 07-6.31 


3.76 


1.36(16) 


0.60(15) 


0.15(18) 


0.63(14) 


1.16(14) 


W3(OH) 


5.09 


1.96(06) 


1.26(08) 


1.07(08) 


1.24(10) 


1.81(09) 


GL490 


3.37 




-0.41(17) 








S231 


4.03 


1.49(16) 


0.96(17) 


0.29(17) 


0.78(15) 


1.37(15) 


S235 


3.72 


1.36(12) 


0.70(13) 


0.31(15) 


0.66(12) 


1.17(10) 


S241 


4.14 


1.49(14) 


0.83(27) 




0.66(26) 


1.39(22) 


S252A 


4.56 


1.41(65) 


0.35(15) 


-0.43(19) 


0.14(15) 


1.19(23) 


S255 


4.20 




0.76(11) 


0.71(18) 






RCW142 


4.82 


2.22(11) 


1.68(12) 


0.97(06) 


1.26(10) 


1.73(09) 


W28A2(1) 


5.45 


1.72(09) 




0.87(05) 


1.64(10) 


1.66(07) 


M8E 


4.26 


1.38(39) 


0.95(08) 


-0.01(13) 


0.64(11) 


1.12(19) 


9.62+0.19 


5.75 


2.34(11) 


2.01(08) 


1.87(20) 


1.75(10) 


2.22(23) 


8.67-0.36 


5.12 






0.84(10) 


1.48(11) 


2.01(29) 


W31(l) 


6.22 




2.57(15) 


2.18(35) 


2.38(11) 




10.60-0.40 


6.05 


2.78(06) 


2.31(06) 


2.42(06) 


2.42(10) 


2.54(08) 


12.42+0.50 


4.27 


1.46(17) 








1.35(21) 


12.89+0.49 


4.66 






0.07(12) 


1.03(12) 




12.21-0.10 


5.98 


3.03(53) 






2.24(16) 


3.16(37) 


W33cont 








3.05(06) 


2.54(11) 




13.87+0.28 


5.29 


2.89(10) 




1.08(19) 


1.15(16) 


2.74(09) 


W33A 


5.19 






0.69(19) 


1.31(11) 




14.33-0.64 


4.33 




1.16(20) 


0.94(17) 


1.11(10) 




19.61-0.23 


5.29 




1.54(19) 


1.04(12) 


1.37(14) 


1.74(59) 


20.08-0.13 


4.65 






0.86(14) 


0.89(11) 




23.95+0.16 


5.36 


2.23(54) 




1.21(22) 


1.42(14) 


2.06(51) 


24.49-0.04 


4.74 






0.42(69) 


0.76(12) 




W42 








1.43(07) 


1.94(13) 




28.86+0.07 


5.53 








1.55(10) 




W43S 


6.24 


2.88(12) 


2.35(12) 


1.19(12) 


2.05(12) 


2.60(09) 


31.41+0.31 


5.32 






1.17(09) 


1.85(11) 




W43 Main3 








1.09(14) 


1.87(18) 




31.44-0.26 


5.78 






1.61(18) 


1.84(11) 




32.05+0.20 


5.02 






1.48(43) 


1.50(11) 




32.80+0.20 








1.93(12) 


2.35(12) 




W44 


5.48 


2.17(21) 


1.90(09) 


1.73(07) 


1.96(11) 


2.26(06) 


S76E 


4.47 






0.83(12) 


1.14(10) 




35.58-0.03 


4.68 


1.40(38) 




0.43(19) 


0.93(15) 


1.41(13) 


35.20-0.74 


4.65 






1.06(18) 


1.27(13) 




W49N 


6.69 


3.44(06) 




3.29(11) 


2.91(11) 


3.22(06) 


W49S 








1.81(51) 


2.44(19) 




43.80-0.13 


4.51 








0.69(11) 




45.07+0.13 




2.42(16) 


2.08(17) 


1.84(18) 


2.04(11) 


2.24(15) 


48.61+0.02 


6.06 


3.12(12) 


2.07(21) 


0.57(22) 


1.73(15) 


2.69(27) 


W51W 


5.93 






1.89(29) 


1.76(23) 
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Table 7 — Continued 



Source 


logL/fl 




L (HOJN3-2j 


L (Ob7-D) 




L (Cb2-1) 




(L©) 


K kms~^pc^ 


K kms pc 


K kms pc 


K kms~^pc'^ 


K kms pc 


VV tJJ-1V± 


6.60 




9 7Q('9n'l 








59.78+0.06 


4.10 


1.53(09) 






0.58(16) 


1.33(14) 


S87 


4.53 


1.74(06) 








1.41(08) 


S88B 


5.03 


1.69(16) 


0.89(17) 


0.47(15) 


0.60(12) 


1.08(10) 


K3-50 


6.32 


2.68(08) 






1.95(15) 


2.32(13) 


ONI 


5.16 


2.26(09) 




1.18(20) 


1.45(14) 


2.14(10) 


ON2S 


5.71 


2.42(14) 


2.21(17) 


2.13(13) 


1.96(12) 


2.33(13) 


ON2N 








1.73(22) 


1.67(13) 




S106 


5.71 


2.26(09) 


1.64(12) 


1.25(44) 


1.13(21) 


2.02(11) 


W75N 


5.29 


2.11(07) 


1.66(06) 


1.50(08) 


1.56(11) 


1.90(11) 


DR21S 


5.70 


0.29(06) 


0.99(05) 


0.81(11) 


1.43(12) 


2.11(10) 


W75(OH) 


4.70 


2.39(09) 


1.64(19) 


1.64(10) 


1.54(11) 


2.24(08) 


BFS 11-B 


3.88 






-0.40(31) 


0.03(39) 




CEPA 


4.92 


0.93(12) 


0.47(10) 






0.91(14) 


S158 


5.29 


2.40(07) 


1.71(12) 




1.63(10) 


2.25(08) 


S158A 






1.67(08) 








NGC7538-IRS9 






1.31(29) 








S157 


4.35 


1.56(14) 


0.84(22) 


0.60(32) 


0.73(13) 


1.25(22) 


IRAS19094 


4.47 


1.86(45) 








1.38(10) 


IRAS20050 


2.47 


0.38(39) 


-0.27(48) 






0.25(70) 


IRAS20106 


3.16 


0.48(34) 








0.23(30) 


IRAS20126 


4.78 


1.74(09) 








1.35(17) 


IRAS20188 


2.48 


-0.15(10) 


-0.76(11) 






-0.24(14) 


IRAS20216 


3.47 


0.65(14) 








0.24(19) 


IRAS20220 


4.78 


2.02(25) 








1.94(81) 


IRAS20333 


3.96 


1.23(08) 


0.54(47) 






0.77(13) 


IRAS21391 


2.43 


0.41(14) 


-0.08(26) 






0.35(11) 


IRAS22172 


3.56 


0.99(20) 








0.70(30) 


IRAS22198 


2.97 


0.63(27) 








0.56(25) 


IRAS23011 


1.95 


0.11(18) 


-0.81(28) 






-0.11(28) 


IRAS23385 


4.73 


2.01(17) 








1.77(17) 


L483 


1.11 




-1.91(89) 








IRAS03282 


0.08 


-0.81(53) 


-2.37(98) 






-0.99(01) 


L1251B 


1.04 


-0.80(17) 








-0.75(35) 


Mean'^ 


4.7 X 10"' 


359(103) 


79(21) 


64(24) 


77(18) 


272(94) 


median"^ 


1.06 X 10^ 


148 


43 


16 


29 


91 



Values in this tabic arc in log scale, except the mean and median at the bottom. Numbers in the parentheses are 
uncertainties of the line luminosity in the log scale, in the unit of 0.01. 

''Converted from Shirley et al. (2003). 

''Statistics on massive dense clumps only, which does not include the sources below the horizontal line. Values and 
uncertainties (in the parentheses) are in the unit of Kkms~^pc^. 
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Table 8: Ratios of Sizes and Luminosities for Different Tracers" 





Hhcni-o 

Rhcnz-2 


RCS2-1 
RcS7-6 


Rhcni-o 

RCS2-1 


RcS5-4 

Rcsj-a 


RhCN3-2 

RcS5-4 


Mean'' 


2.97±0.19 


4.12±0.45 


1.11±0.03 


1.42±0.09 


1.11±0.04 


Median 


2.71 


3.36 


1.07 


1.31 


1.05 


Sdv'^ 


1.06 


2.40 


0.23 


0.64 


0.21 




Tl 

^HCNl-0 


Tl 

'-'CS2-1 


Tl 

^CS2-1 


Tl 


Tl 

^HCN3-2 




tI 

^HCN3-2 


Tl 


Tl 


Tl 

^CS7-6 


Tl 

^CS5-4 


mean'' 


4.3±0.4 


9.4±2.1 


4.0±0.4 


2.6±0.4 


1.6±0.1 


median 


3.8 


5.9 


3.4 


1.7 


1.5 


Sdv^ 


2.5 


11.0 


2.3 


2.6 


0.7 



"Sizes of CS 5-4 maps come from Shirley et al. (2003). 

'The uncertainty is the sigma of the mean 
'^The standard deviation of the ratios 



Table 9: Properties of Different Transitions 



Transition 




'-'up 




Mean L' 


Mean L;„,- 


Mean R 


Mean Kgub'^ 




(lOOK) 


(K) 


(lOOK) 


Kkm/spc^ 


Kkm/spc^ 


(pc) 


(pc) 


CS 7-6 


2.0x10'^ 


66 


2.6x10^ 


64 


42 


0.28 


0.24 


CS 5-4 '' 


6.9x10^ 


35 


6.0x10"^ 


77 


68 


0.37 


0.31 


HCN 3-2 


6.8x10'^ 


26 


3.6x10^ 


79 


101 


0.32 


0.35 


CS 2-1 


3.9x10^ 


7.1 


4.1x10=^ 


272 




1.03 




HCN 1-0 


4.5x10^ 


4.3 


5.1x10^ 


359 




1.13 





"Prom Evans (1999). 
'From Shirley et al.( 2003). 

■^Statistics on a subsample of the sources only, which have been mapped by all the 3 high-J transitions CS 7-6, CS 5-4, 
and HCN 3-2 



Table 10: Statistics of linewidth for different tracers 



Tansitions 


Ay(CS7-6) 


AF(CS5-4) " 


Ay(CS2-l) 


Ay(HCN 1-0) 


AF(C-^'iS5-4) 




(km s"""^) 


(km s~^) 


(km s~^) 


(F=2-l)(km s"^) 


(km s"""^) 


Mean 


5.9±2.3 


5.6±2.2 


4.8±1.8 


3.9±1.0 


4.7±1.9 


]\l(xliaii 


5.7 


5.1 


1.8 


3.7 


1.2 



"Prom Shirley et al.( 2003). 



Table 11. The Virial Mass, Surface Density and Average Volume Density of Dense Clumps ^ 



SOURCE 


^^lvlQ 


HON 3^2 

s 






HON 1 

( EF CTn~'^) 

V& ^^^^ J 






CS 7 -J- 6 
( er rm— 


loci Tl 1 


(Mr,) 


CS 2 ^ 1 

(^ CTn~'^) 
V& ^^^^ } 




G121. 30+0.66 


2.17(12) 


0.69 


5.54 


2.35(10) 


0.45 


5.19 


1.79(16) 


1.65 


6.31 


2.44(11) 


0.37 


5.01 


G123.07-6.31 


2.68(16) 


2.58 


6.15 


3.24(16) 


0.71 


5.03 


2.52(17) 


3.75 


6.48 


3.21(15) 


0.76 


5.09 


W3(OH) 


2.99(05) 


1.40 


5.60 


3.43(03) 


0.51 


4.72 


2.74(05) 


2.54 


6.11 


3.39(05) 


0.56 


4.80 


S231 


2.24(10) 


0.33 


5.03 


2.60(09) 


0.14 


4.30 


2.02(10) 


0.55 


5.47 


2.61(09) 


0.14 


4.29 


S235 


2.08(09) 


0.38 


5.21 


2.49(09) 


0.15 


4.38 


1.90(10) 


0.58 


5.57 


2.42(08) 


0.18 


4.53 


S255 


2.25(08) 


0.71 


5.53 








2.15(10) 


0.91 


5.74 








RCW142 


2.97(06) 


1.60 


5.70 


3.45(06) 


0.53 


4.74 


2.73(03) 


2.73 


6.16 


3.23(05) 


0.87 


5.17 


W28A2(1) 








3.28(04) 


0.95 


5.20 


2.70(01) 


3.66 


6.37 


3.20(04) 


1.15 


5.36 


M8E 


2.07(05) 


0.32 


5.10 


2.55(18) 


0.10 


4.13 


1.78(08) 


0.61 


5.67 


2.42(10) 


0.14 


4.40 


G9.62+0.19 


3.49(04) 


1.01 


5.14 


3.89(06) 


0.40 


4.34 


3.39(10) 


1.26 


5.33 


3.86(12) 


0.43 


4.40 


G8.67-0.36 


3.17(06) 


0.68 


5.04 








2.84(06) 


1.45 


5.70 








W(31) 


3.89(09) 


1.05 


4.97 








3.96(18) 


0.89 


4.82 








GlO.60-0.40 


3.41(02) 


1.19 


5.28 


3.85(03) 


0.44 


4.42 


3.46(03) 


1.06 


5.18 


3.80(04) 


0.49 


4.51 


G12.21-0.10 








4.63(24) 


0.09 


2.98 






7.42 


4.69(16) 


0.08 


2.85 


G12.89+0.49 














2.32(07) 


1.49 


5.98 








W33cont 














3.36(03) 


0.45 


4.68 








G13. 87+0.28 








3.22(08) 


0.04 


3.10 


2.41(10) 


0.23 


4.71 


3.09(06) 


0.05 


3.35 


G14.33-0.64 


2.18(11) 


0.57 


5.42 








2.25(10) 


0.48 


5.27 








G23.95+0.16 








3.12(23) 


0.04 


3.18 


2.51(13) 


0.15 


4.40 


3.09(22) 


0.04 


3.24 


W43S 


3.12(07) 


0.29 


4.50 


3.64(06) 


0.09 


3.46 


2.76(07) 


0.64 


5.21 


3.53(05) 


0.11 


3.68 


G31.41+0.31 














3.11(05) 


1.39 


5.54 








G31. 44-0.26 














3.09(11) 


0.26 


4.45 








G32.05+0.06 














3.47(21) 


0.70 


4.91 








W44 


2.97(09) 


0.42 


4.82 


3.49(13) 


0.13 


3.78 


2.79(08) 


0.63 


5.18 


3.29(08) 


0.20 


4.18 


G35.58-0.03 








3.51(21) 


0.54 


4.72 


2.86(14) 


2.41 


6.02 


3.46(11) 


0.61 


4.82 


G35.20-0.74 














2.81(10) 


1.13 


5.55 








W49S 














3.73(24) 


0.93 


4.96 








G45.07+0.13 


3.52(10) 


0.71 


4.90 


3.90(09) 


0.30 


4.13 


3.39(10) 


0.95 


5.15 


3.66(09) 


0.52 


4.62 


G48. 61+0.02 


2.73(13) 


0.08 


3.88 


3.47(07) 


0.02 


2.41 


2.31(14) 


0.22 


4.72 


3.40(14) 


0.02 


2.54 


W51W 














2.96(15) 


0.22 


4.40 








W51M 


3.98(10) 


1.00 


4.89 


4.51(17) 


0.29 


3.82 


3.75(03) 


1.68 


5.34 


4.44(03) 


0.34 


3.96 


G59.78+0.06 








2.53(10) 


0.11 


4.14 








2.54(11) 


0.10 


4.13 


S87 








2.60(06) 


0.12 


4.18 








2.57(06) 


0.13 


4.24 



Table 11 — Continued 



SOURCE 




HON 3^2 






HON 1^0 






CS 7 ^ 6 






OS 2 -5> 1 






log(M„i^) 


E 


log(n) 




S 


log(n) 




E 


log(n) 


log(M„i^) 


E 


log(n) 




(M©) 


(g cm-2) 




(M©) 


(g cm-2) 




(M©) 


(g cm-2) 




(M©) 


(g cm-2) 




S88B 


2.15(10) 


0.33 


5.07 


2.69(09) 


0.10 


4.00 


2.06(09) 


0.41 


5.27 


2.43(06) 


0.17 


4.52 


K3-50 








4.12(07) 


0.38 


4.18 








4.09(09) 


0.40 


4.23 


ONI 








3.54(06) 


0.14 


3.80 


2.81(12) 


0.73 


5.26 


3.49(07) 


0.15 


3.91 


ON2S 


3.07(09) 


0.23 


4.37 


3.34(08) 


0.12 


3.84 


2.97(07) 


0.29 


4.58 


3.34(07) 


0.12 


3.84 


ON2N 














2.86(11) 


0.40 


4.85 








W75N 


2.85(07) 


0.66 


5.18 


3.22(03) 


0.28 


4.44 


2.72(04) 


0.88 


5.43 


3.18(06) 


0.31 


4.53 


DR21S 


2.93(03) 


1.03 


5.43 


3.43(03) 


0.33 


4.43 


2.96(06) 


0.97 


5.38 


3.45(05) 


0.31 


4.39 


W750H 


3.22(10) 


0.79 


5.12 


3.59(05) 


0.34 


4.37 


2.97(06) 


1.40 


5.61 


3.54(04) 


0.37 


4.46 


CepA 


2.24(06) 


1.56 


6.05 


2.67(07) 


0.58 


5.19 








2.73(07) 


0.50 


5.07 


S158 


3.11(06) 


0.40 


4.72 




















S158A 


2.87(06) 


0.35 


4.76 




















IRAS20050+2720 


2.06(29) 


1.32 


6.03 


2.53(26) 


0.45 


5.08 








2.52(51) 


0.46 


5.11 


IRAS20106+3545 








2.17(24) 


0.17 


4.64 








2.22(23) 


0.15 


4.55 1 


IRAS20126+4104 








2.83(15) 


0.20 


4.40 








2.84(17) 


0.19 


4.3^ 


IRAS20188+3928 


1.75(17) 


5.22 


7.08 


2.19(17) 


1.89 


6.20 








2.22(17) 


1.76 


6.1SP 


IRAS20216+4107 


1.49(34) 


0.96 


6.10 


2.24(20) 


0.17 


4.62 








2.14(21) 


0.22 


1 

4.81 ' 


IRAS20220+3728 








3.43(20) 


0.09 


3.61 








3.44(30) 


0.09 


3.60 


IRAS20333+4102 


2.14(25) 


0.12 


4.43 


2.33(14) 


0.08 


4.06 








2.30(15) 


0.08 


4.11 


IRAS21391+5802 


1.92(20) 


1.42 


6.14 


2.27(16) 


0.63 


5.44 








2.25(15) 


0.66 


5.47 


IRAS22172+5549 








2.56(21) 


0.19 


4.51 








2.60(24) 


0.17 


4.42 


IRAS22198+6336 








2.29(21) 


0.13 


4.42 








2.32(20) 


0.13 


4.37 


IRAS23011+6126 


1.29(25) 


1.40 


6.45 


2.05(21) 


0.24 


4.93 








2.03(24) 


0.25 


4.97 


IRAS23385+6053 








3.41(18) 


0.13 


3.86 








3.38(18) 


0.14 


3.91 


IRAS03282+3035 


-0.09(32) 


0.09 


5.34 


0.52(21) 


0.02 


4.13 








0.48(24) 


0.02 


4.19 


L483 


-0.06(27) 


0.07 


5.44 




















L1251B 








1.00(09) 


0.1.') 


'i.r.i 








i.iUio) 


0.11 


l.«5 


mean'' 


1555 


0.78 


2.5x10^ 


5300 


0.29 


3.2x10"* 


1300 


1.09 


1.7xlO«* 


5000 


0.33 


4.2x10'' 


sigma ^ 


445 


0.55 


0.7x105 


1760 


0.23 


0.8 xlO^ 


310 


0.91 


l.OxlO^ 


1870 


0.27 


1.0 xlO* 


median'' 


890 


0.68 


1.3x10^ 


2700 


0.28 


1.6x10-* 


650 


0.89 


3.7x10^ 


1950 


0.31 


2.6x10-* 



^Thc virial mass, surface density and volume density are calculated for mass within the FWHM contour of each transitions. Values for virial mass and volume density are 
in log scale, except the mean and median at bottom. Numbers in the parentheses are uncertainties in the log scale, in the unit of 0.01. 

''Statistics on massive dense clumps only, which does not include the sources below the horizontal line in the table. 



"^The c of the virial mass is the c of the mean. 



I 

CO 

I 



-54- 



Table 12: Statistics of Line Ratio log S'^^ by Evolutionary Stages " 





No HII region 


UC HII 


Comp. HII 


HII 


mean'' 


-0.45 ±0.08 


-0.21±0.08 


-0.40±0.15 


-0.21±0.08 


median 


-0.47 


-0.20 


-0.28 


-0.21 


Std. dev.^ 


0.34 


0.37 


0.41 


0.25 



"uminosities of CS 5-4 come from Shirley et al. (2003). 

'The uncertainty is the sigma of the mean 
'^The standard deviation of the ratios 
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Table 13. Evolutionary Analysis of Massive Clumps 



Source 


HII? 


[o„ — j±m ^ 


b,c in„i'(«CAr3-2) b 


l„„i'(CS7-6) b,c 


Outflow? d Infall? " 


'"8M„i^(«CS5-4) 

(M©j 


S L' (HCNl—0) 


°L'(CSB-4) 


G121. 30+0.66 




0.34 


-0.44 


-0.47 


Y Y 


S231 




1.77 


-0.53 


-0.49 


Y 


S241 




1.99 


-0.66 




Y 


12.89+0.49 




2.13 




-0.96 




W33A 




2.53 




-0.62 




24.49-0.04 




2.26 




-0.34 




28.86+0.07 




2.68 








W43 Main3 








-0.78 




32.05+0.20 




1.85 




-0.02 


Y 


W75(OH) 




1.60 


-0.66 


0.10 


Y 


BFS 11-B 




1.84 




-0.43 


Y 


Mean* 




2.02(0.12) 


-0 62^0 04") 


-0 45(0 1 1 ■) 




Median 




1.99 


-0.66 


-0.47 






UCHII 


2.08 


-0.70 


-0.17 


Y 


S255 


UCHII 








Y 


Rr,W142 


UCHII 


2.25 


-0.54 


-0.29 


Y 


W28A2('1'I 

VV £j\DIi^^ I J- J 


UCHII 


2.80 




-0.77 


Y 


M8E 


UCHII 


2.26 


-0.43 


-0.65 


Y 


q fiO-l-O 1 Q 


UCHII 


2.40 


-0.24 


0.12 


Y 


8 67-0 36 


UCHII 


2.19 




-0.64 


Y 


W31 n ~i 


UCHII 


2.36 




-0.20 


Y 


1 2 42-4-0 "iO 


UCHII 








Y 


VV tJOL/'JlHj 


UCHII 






0.51 




1 4 33-0 64 


UCHII 


2.13 




-0.17 




20 08-0 1 3 


UCHII 


1.44 




-0.03 




W42 


UCHII 






-0.51 




31.41+0.31 


UCHII 


2.30 




-0.68 


Y 


31.44-0.26 


UCHII 


2.73 




-0.23 




35.58-0.03 


UCHII 


1.57 




-0.50 




W49N 


UCHII 


2.53 




0.38 


Y 


W49S 


UCHII 






-0.63 




43.80-0.13 


UCHII 


2.08 








45.07+0.13 


UCHII 




-0.34 


-0.20 


Y 


59.78+0.06 


UCHII 


2.62 






Y Y 


S87 


UCHII 








Y 


S88B 


UCHII 


2.83 


-0.80 


-0.13 


Y 


ONI 


UCHII 


2.04 




-0.27 




S106 


UCHII 


2.85 


-0.62 


0.12 




W75N 


UCHII 


2.44 


-0.43 


-0.06 


Y 


DR21S 


UCHII 


2.70 


-0.30 


0.38 


Y Y 


CEPA 


UCHII 




-0.46 




Y Y 


S158 


UCHII 


2.33 


-0.69 




Y Y 


Mean' 




2.32(0.08) 


-0.48(0.06) 


-0.21(0.08) 




Median 




2.35 


-0.45 


-0.20 




10.60-0.40 


CHII 


2.61 


-0.38 




Y 



Table 13— Continued 



-feziy ' l°g^S^ Outflow? d Infall? « 

(Mo) 



12.21-0.10 


CHII 


2.30 










19.61-0.23 


CHII 


2.19 




-0.33 






W43S 


CHII 


3.21 


-0.53 


-0.86 




Y 


32.80-I-0.20AB 


CHII 






-0.42 






W44 


CHII 


2.33 


-0.18 


-0.23 




Y 


48.6H-0.02 


CHII 


3.42 


-1.05 


-1.16 






M51M 


CHII 


2.83 


-0.61 


-0.11 




Y 


K3-50 


CHII 


2.53 






Y 




ON2N 


CHII 






0.06 


Y 




S157 


CHII 


2.05 


-0.72 


-0.13 




Y 


Mean' 




2.61(0.18) 


-0.62(0.14) 


-0.40(0.15) 






Median 




2.43 


-0.61 


-0.28 






G123.07-6.3 


HII 


1.06 


-0.76 


-0.48 


Y 


Y 


S235 


HII 


1.72 


-0.66 


-0.35 






S252A 


HII 


2.41 


-1.06 


-0.57 






13.87-1-0.28 


HII 


2.80 




-0.07 






23.95+0.16 


HII 


2.73 




-0.21 




Y 


S76E 


HII 


2.09 




-0.31 


Y 




35.20-0.74 


HII 


1.14 




-0.21 






W51W 


HII 


2.89 




0.13 






ON2S 


HII 


2.62 


-0.21 


0.17 






Mean* 




2.30(0.22) 


-0.64(0.25) 


-0.21(0.08) 






Median 




2.52 


-0.66 


-0.21 







f 

The mean in log scale. The uncertainty is the sigma of the mean. 

''The classification of evolutionary status of HII regions are judged from the size of the HII region from 2cm 
survey. Refer to Shirley at al. 2003 and references in there for details. 

''Values are in log scale. 

"^CS 5-4 data from Shirley et al. 2003. 

''Refer to outflow catalogue, Wu et al. (2004). 

''Judged from HCN 3-2 proflle for Sv < -0.25. 



Tabic 14: Statistics of Ratio ^"'^^^ 

molcc ulc 

log ^Vir(RcS2-ll log ^Viri^CSr-fil Ipg ^VirjRHCNl^ log ^^^"^"^^ 
^CS2-1 ^CS7-6 ^HCNl-0 ^ 



mean" 


1.34 ±0.07 


1.57±0.07 


1.30±0.10 


1.35±0.07 


median 


1.30 


1.59 


1.20 


1.29 


Std. dev.'' 


0.39 


0.43 


0.54 


0.37 



"The uncertainty is the sigma of the mean 
''The standard deviation of the ratios 
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Table 15: Statistics of Properties of Massive Dense Clumps " 





^vir{HCNl-0) 

(Mo) 


(Mo) 


Mmr(C52-l) 

(Mo) 


M„ir(CS7-6) 

(Mo) 


L/i?, 
(Lo) 


Mean^ 


5300±1760 


1555±445 


5000±1870 


1300±310 


470000± 120000 


Median 


2700 


890 


1950 


650 


106000 


Sdv'' 


9600 


2310 


10300 


1930 


920000 










^IR 






^vir(HCNl-O) 

(Lo /Mo) 


^vir(,HCN3-2) 

(Lo /Mo) 


^vir(CS2-l) 

(Lo /Mo) 


^vir(CS7-6) 

(Lo /Mo) 




mean'' 


117±20 


371±102 


139±26 


592±195 




median 


98 


182 


129 


229 




Sdv'^ 


105 


487 


134 


1085 





"Statistics on massive dense clumps only, which does not include the sources below the horizontal line in Table 1. 
'The uncertainty is the sigma of the mean. 
■^The standard deviation. 



Table 16: Statistics of Properties of Massive Dense Clumps with Lm > 10^^^ Lo " 





^vir{HC'Nl-0) 

(Mo) 


^vir{HCN3-2) 

(Mo) 


Mto,.(C52-1) 

(Mo) 


M„jr.(C57-6) 

(Mo) 


L//? 
(Lo) 


Mean'' 


6650±2350 


2070±630 


6300ib2500 


1480±420 


630000± 160000 


Median 


2950 


1170 


2500 


710 


200000 


Sdv'^ 


10800 


2600 


11500 


2030 


1020000 




^IR 


J^IR 


J^IR 


i^IR 






^vir(HCNl-a) 

(L„. /M,J 


Mr,r(JICN:i-2) 

(L„. /M,J 


M,jir(CS2-l) 


(L,. /M,J 




mean'' 


137±23 


473±130 


169±30 


726±246 




median 


120 


324 


155 


439 




Sdv^ 


106 


532 


137 


1204 





"Statistics on massive dense clumps only, which does not include the sources below the horizontal line in Table 1. 
'The uncertainty is the sigma of the mean. 
■^The standard deviation. 




Fig. 1. — HCN 1-0 contour maps of massive clumps. The lowest contour level and increasing step 
of contours are indicated in the plot. The beam size is shown at the lower left of each map. 
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Fig. 2. — HCN 1-0 contour maps continue. 
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Fig. 3. — HCN 1-0 contour maps continue... 




Fig. 4. — HCN 1-0 contour maps continue. 
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Fig. 5. — HCN 1-0 contour maps continue. 
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Fig. 6. — HCN 1-0 contour maps continue... 
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Fig. 7. — HCN 1-0 contour maps continue. 




Fig. 8. — HCN 3-2 contour maps of massive clumps. The lowest contour level and increasing step 
of contours are indicated in the plot. The beam size is shown at the lower left of each map. 




Fig. 9. — HCN 3-2 contour maps continue. 




Fig. 10. — HCN 3-2 contour maps continue. 




Fig. 11. — HCN 3-2 contour maps continue. 




Fig. 12. — CS 2-1 contour maps of massive clumps. The lowest contour level and increasing step 
of contours are indicated in the plot. The beam size is shown at the lower left of each map. 




Fig. 13. — CS 2-1 contour maps continue. 
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Fig. 14. — CS 2-1 contour maps continue. 




Fig. 15.— 



CS 2-1 contour maps continue. 




Fig. 16. — CS 2-1 contour maps continue. 
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Fig. 17. — CS 2-1 contour maps continue. 
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Fig. 18. — CS 2-1 contour maps continue. 




Fig. 19. — CS 7-6 contour maps of massive clumps. The lowest contour level and increasing step 
of contours are indicated in the plot. The beam size is shown at the lower left of each map. 




Fig. 20. — CS 7-6 contour maps continue. 




Fig. 21. — CS 7-6 contour maps continue. 




Fig. 22. — CS 7-6 contour maps continue. 
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Fig. 23. — CS 7-6 contour maps continue. 




Fig. 24. — The plots to show the effect when we smooth the HCN 3-2 maps (left panels) of four 
clumps with more than one peaks, to the resolution of HCN 1-0 data (middle panels), comparing 
to the real HCN 1-0 maps (right panels). 
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Fig. 25. — The one-dimensional radial profiles of some dense clumps. The intensities have been 
normalized to the peak values, and the angular radii (convolved with beam sizes) are in units of 
the FWHM radius of each source. 
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Fig. 26. — Histograms of the deconvolved angular size of massive dense clumps in different transi- 
tions. The vertical dashed line in each panel indicates the beam size of each transitions. 
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Fig. 27. — Histograms of the deconvolved linear size of massive dense clumps in different transitions. 
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Fig. 28. — Histograms of the infrared and line luminosity of massive dense clumps in different 
transitions. 
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Fig. 29. — The inverse AV-RpwHM relation between different tracers for massive dense clumps. 




Fig. 30. — Correlations between the infrared luminosity and the virial mass measured by various 
dense gas tracers. The fitted lines show the linear least squares fit for clumps with Ljr > lO^^^ Lq. 
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Fig. 31. — Plots of the surface density vs. the ratio of infrared luminosity to virial mass, infrared 
luminosity, line luminosity, and virial mass derived from HCN 3-2 maps. 
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Fig. 32. — Plots of the surface density vs. the ratio of infrared luminosity to virial mass, infrared 
luminosity, line luminosity, and virial mass derived from HCN 1-0 maps. 
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Fig. 33. — Plots of the surface density vs. the ratio of infrared luminosity to virial mass, infrared 
luminosity, line luminosity, and virial mass derived from CS 7-6 maps. 
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Fig. 34. — Plots of the surface density vs. the ratio of infrared luminosity to virial mass, infrared 
luminosity, line luminosity, and virial mass derived from CS 2-1 maps. 



-92- 



C7> 
O 



O 



I I I I I I I I I I I I I I I I I I I I I I I I 



V V ' 




0! 



o 



a: 

j; 
o 



4 6 
log(<n>cs7-6) 



- • • 



I I I I I I I I I I I I I I I I I I I I I 



3 4 5 

log(<n>HCNi_o) 




4 6 
log(<n>Hcio-2) 



Fig. 35. — Correlation between infrared luminosity and the mean volume density (n) (n in the text) 
within the FWHM contour of CS 2-1, CS 7-6, HCN 1-0, and HCN 3-2 maps. The dashed line has 
the slope predicted by Krumholz & McKee (2008) and Narayanan et al. (2008) and an arbitrary 
offset. The clumps with Lm above and below the threshold luminosity lO^'^L© are plot with solid 
circles and hollow stars, respectively. 
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Fig. 36. — Correlation between line luminosity and the mean volume density (n) {n in the text) 
within the FWHM contour of CS 2-1, CS 7-6, HCN 1-0, and HCN 3-2 maps. 
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Fig. 37. — Liji vesus volume density derived from LVG model (Plume et al. (1997). 
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Fig. 38. — The star formation efficiency indicated by the infrared luminosity per unit of dense gas 
mass vs. virial mass. 
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Fig. 39. — Luminosity ratios between different J transitions versus infrared luminosity. 
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Fig. 40. — Luminosity ratios between different J transitions versus bolometric temperature. 
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Fig. 41.— 



Luminosity ratios between different J transitions versus star formation efficiency. 



-99- 



; • Hii 


1 




- X CHII 






K UCHI 






- QNoRC 




n u 


o 


• 
• 


JO 

K O 




1 


1 < < < < 1 < < < 



12 3 4 

log(L|R) / log(MM^i,(^35_4)) 



Fig. 42. — The distribution of clumps at different evolutionary states in the line luminosity ratio 
versus star formation efficiency plot. 
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Fig. 43. — Correlation between the line luminosity of different transitions and the virial mass within 
FWHM of contour maps. The dashed-line shows the linear least squares fit with uncertainties; the 
solid line indicates the robust fit. Data of CS 5-4 comes from Shirley et al. (2003). 
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Fig. 44.— LjR-L[^^^^i correlation for HCN 1-0, HCN 3-2, CS 2-1, CS 5-4 and CS 7-6 Galactic clumps. 
Data of CS 5-4 are taken from Shirley et al. 2003. The squares in the upper left panel are HCN 
1-0 observations of gaIaxies(Gao &; Solomon 2004a) with 3 high-z detections ( Solomon et al. 2003, 
Vanden Bout et al. 2004, and Carilli et al. 2005). The squares in the upper right panel show the 
HCN 3-2 observations of galaxies (Bussmann et al. 2008). The solid line in the first panel shows 
the best fit Liji-L'fj^j^^_Q correlation from galaxies, and is shown in other plots to indicate the 
correlation shifts between HCN 1-0 and other tracers. The dashed line shows where L7i^=lO^'^L0 
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Fig. 45. — Correlations between the distant independent ratio hj^/L'^^^ vs. Ljr for different 
tracers. The squares in figures are HCN 1-0 observations of galaxies (Gao &: Solomon 2004a) in 
the first panel, and HCN 3-2 observations of galaxies (Bussmann et al. 2008) in the second panel. 
The horizontal dashed line in the top plot indicates the averaged ^ir/ L'jjq^-^^_q ratio for galaxies; 
the vertical dashed line in the top plot shows the cutoff at 'Lir=1Q'^-^'Lsun- These two lines are 
also shown in other plots to indicate the relative shifts of Ljr/L'^^i between HCN 1-0 and other 
tracers. 
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Fig. 46. — Log-log scale plots of the surface density vs. distance for massive clumps for different 
transitions. The lines in the plots have fixed slopes of —1 with different offsets. 
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Fig. 47. — Log-log scale plots of the areraged volume density vs. distance for massive clmiips with 
different transitions. The Unes in the plots have fixed slopes of —2 with different offsets. 
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Fig. 48. — Plots of the line ratios vs. distance for massive clumps. No obvious trends are seen. 



